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as well as some of the novel applications of graphene. In particular, this review cov- 
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properties, magneto-transport, integer and fractional quantum Hall effects, mechanical 
properties, transistors, optoelectronics, graphene-based sensors, and biosensors. This 
approach attempts to highlight both the means by which the current understanding of 
graphene has come about and some tools for future contributions. 
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Graphene is a single two-dimensional layer of carbon 
atoms bound in a hexagonal lattice structure. It has been 
extensively studied in the last several years even though 



it was only isolated for the first time in 2004 (Novoselov 



et al. 2004). Andre Geim and Konstantin Novoselov 



won the 2010 Nobel Prize in Physics for their ground- 
breaking work on graphene. The fast uptake of interest 
in graphene is due primarily to a number of exceptional 
properties that it has been found to possess. 

There have been several reviews discussing the topic 
of graphene in recent years. Many are theoretically ori- 
ented, with Castro Neto et a/.'s review of the electronic 



properties as a prominent example (Castro Neto et al 



2009 ) and a more focused review of the electronic trans- 
port properties ( [Pas Sarma et a/. , 2011). Experimental 
reviews, to name only a few, include detailed discussions 



of synthesis (Choi et a/., 2010) and Raman character- 
ization methods (Ni et al. 2008a), of transport mech- 
anisms (Avouris, 2010 Giannazzo et al. 2011), of rel- 



evant applications of graphene such as transistors and 
the related bandgap engineering fSchwierzl|2010D, and of 



graphene optoelectronic technologies (Bonaccorso et al. 



2010). We feel, however, that the literature is lacking a 



comprehensive overview of all major recent experimen- 
tal results related to graphene and its applications. It 
is with the intent to produce such a document that we 
wrote this review. We gathered a great number of re- 
sults from what we believe to be the most relevant fields 
in current graphene research in order to give a starting 
point to readers interested in expanding their knowledge 
on the topic. The review should be particularly well- 
suited to graduate students who desire an introduction 
to the study of graphene that will provide them with 
many references for further reading. 



We have attempted to deliver an up-to-date account 
of most topics. For example, we present recent results 
on the fractional quantum Hall effect and some of the 
newest developments and device details in optoelectron- 
ics. In addition, we include a summary of the work that 
has been done in the field of graphene biosensors. As 
a practical tool, we also give comparative analyses of 
graphene substrate properties, of available bandgap en- 
gineering techniques and of photovoltaic devices in order 
to provide the researcher with useful and summarized 
laboratory references. 

This review is structured as follows: We start by 
reviewing the electronic band structure and its associ- 
ated properties in Sec. [TT] before introducing the vibra- 
tional properties, including the phonon dispersion, in 



Sec. |III[ We then detail the various synthesis methods for 
graphene monolayers and their corresponding character- 
ization. The main properties are then discussed, starting 
with the electric field effect, which originally spurred the 
intense activity in graphene research. This is then fol- 
lowed by a review of the magneto-transport properties, 
which includes the quantum Hall effect and recent results 
on the fractional quantum Hall effect. We conclude the 
review of the main properties with a discussion of the me- 
chanical properties of graphene, before turning to some 
important applications. The transistor applications are 
discussed first, with particular emphasis on band gap en- 
gineering, before reviewing optoelectronic applications. 
We finish the applications section with a discussion of 
graphene-based sensors and biosensors. 



II. ELECTRONIC STRUCTURE 

Graphene has a remarkable band structure thanks to 
its crystal structure. Carbon atoms form a hexagonal 
lattice on a two-dimensional plane. Each carbon atom 
is about a = 1.42 A from its three neighbors, with each 
of which it shares one a bond. The fourth bond is a 
TT-bond, which is oriented in the z-direction (out of the 
plane). One can visualize the tt orbital as a pair of sym- 
metric lobes oriented along the z-axis and centered on the 
nucleus. Each atom has one of these 7r-bonds, which are 
then hybridized together to form what are referred to as 
the TT-band and 7r*-bands. These bands are responsible 
for most of the peculiar electronic properties of graphene. 

The hexagonal lattice of graphene can be regarded as 
two interleaving triangular lattices. This is illustrated in 
Fig. [l] This perspective was successfully used as far back 
as 1947 when Wallace calculated the band structure for a 
single graphite layer using a tight-binding approximation 



(Wallace, 1947). 



Band structure is most often studied from a standpoint 
of the relationship between the energy and momentum 
of electrons within a given material. Since graphene con- 
strains the motion of electrons to two dimensions, our 
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ley degeneracy of = 2. The K and K' points are the 
primary points of interest when studying the electronic 
properties of graphene. This is noteworthy in comparison 
to traditional semiconductors where the primary point of 
interest is generally F, where momentum is zero. 

Following the original work (Wallace, 1947), the tight 
binding Hamiltonian, written in a basis of sublattice po- 
sitions A and B, is given by: 




FIG. 1 Triangular sublattices of graphene. Each atom 
in one sublattice (A) has 3 nearest neighbors in sublattice (B) 
and vice-versa. 




FIG. 2 First Brillouin zone and band structure of 
graphene. The vertical axis is energy, while the horizon- 
tal axes are momentum space on the graphene lattice. The 
first Brillouin zone of graphene is illustrated in the horizon- 
tal plane and labeled with some points of interest. K and K' 
are the two non- equivalent corners of the zone, and M is the 
midpoint between adjacent K and K' points. F is the zone 
center. K and K' are also known as the Dirac points. The 
Dirac points are the transition between the valence band and 
the conduction band. 



CA te 
c.c. 



ik-a 



te 



(1) 



where = = are the on-site energies of the 
carbon atoms on sites A and B, t o:^ 2.7 eV is the next 
nearest hopping element (between A and B sites), al = 
a(l,0), 02 = a( — 1/2, sin[7r/3]), ds = a( — 1/2, — sin[7r/3]) 
are the positions of the three nearest neighbors, and c.c. 
the complex conjugate of the off-diagonal matrix element. 
The eigenvalues of this tight binding Hamiltonian are 
shown in Fig. [2] as a function of ^ = {kx^ ky). 

A. Massless Dirac fermions 

Looking closely at the region near one of the Dirac 
points (K or K') in Fig. |2j the cone-like linear disper- 
sion relation is evident. The Fermi energy for neutral (or 
ideal) graphene is at the Dirac energy, which is the energy 
of the Dirac point. In graphene devices, the Fermi energy 
can be significantly different from the Dirac energy. 

Electrons within about 1 eV of the Dirac energy have 
a linear dispersion relation. The linear dispersion re- 
gion is well-described by the Dirac equation for massless 
fermions. That is, the effective mass of the charge car- 
riers in this region is zero. The dispersion relation near 
the K points is generally expressed as follows: 



E±(k) ^ ±nvF\k-K\ 



(2) 



which corresponds to the spectrum of the Dirac- 
like Hamiltonian for low-energy massless Dirac fermions 
(again in the sublattice basis {A,B}): 



momentum space is also constrained to two dimensions. 
A plot of the energy versus momentum dispersion rela- 
tion for graphene can be found in Fig. |2] using the tight 
binding approximation discussed below. 

Graphene is a zero-gap semiconductor because the con- 
duction and valence bands meet at the Dirac points (see 
Fig. |2|. The Dirac points are locations in momentum 
space, on the edge of the Brillouin zone. There are two 
sets of three Dirac points. Each set is not equivalent 
with the other set of three. The two sets are labeled K 
and K'. The two sets of Dirac points give graphene a val- 



^^ = ^^Fl^k, + ik^ Q'^')=nv^a-k (3) 

This Dirac Hamiltonian is simply the tight binding 
Hamiltonian from Eq. ([T]) expanded close to K with 
hvF = 3ta/2. Close to K' the Hamiltonian becomes 
Hk' = fivp^* • k^ where a = {(Jx^cFy) is the 2D vector 
of the Pauli matrices (and * denotes the complex con- 
jugate), k is the wavevector, and the Fermi velocity is 
vp ~ 10^ m/s, or l/300th the speed of light in vacuum. 
Charge carriers in graphene behave like relativistic par- 
ticles with an effective speed of light given by the Fermi 
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velocity. This behavior is one of the most intriguing as- 
pects about graphene, and is responsible for much of the 
research attention that graphene has received. 



B. Chirality 

Transport in graphene exhibits a novel chirality which 
we will now briefly describe. Each graphene sublattice 
can be regarded as being responsible for one branch of 
the dispersion. These dispersion branches interact very 
weakly with one another. 

This chiral effect indicates the existence of a pseu- 
dospin quantum number for the charge carriers. This 
quantum number is analogous to spin but is completely 
independent of the 'real' spin. The pseudospin lets us 
differentiate between contributions from each of the sub- 
lattices. This independence is called chirality because 
of the inability to transform one type of dispersion into 
another (Das Sarma a/., 2011). A typical example of 



pending on the angle of incidence to the barrier (Kat 



snelson et al , 2006). 



chirality is that you cannot transform a right hand into a 
left hand with only translations, scalings, and rotations. 
The chirality of graphene can also be understood in terms 
of the Pauli matrix contributions in the Dirac-like Hamil- 
tonian described in the previous section. 



C. Klein paradox 

A peculiar property of the Dirac Hamiltonian is that 
charge carriers cannot be confined by electrostatic po- 
tentials. In traditional semiconductors, if an electron 
strikes an electrostatic barrier that has a height above the 
electron's kinetic energy, the electron wavefunction will 
become evanescent within the barrier and exponentially 
decay with distance into the barrier. This means that 
the taller and wider a barrier is, the more the electron 
wavefunction will decay before reaching the other side. 
Thus, the taller and wider the barrier is, the lower the 
probability of the electron quantum tunneling through 
the barrier. 

However, if the particles are governed by the Dirac 
equation, their transmission probability actually in- 
creases with increasing barrier height. A Dirac electron 
that hits a tall barrier will turn into a hole, and propagate 
through the barrier until it reaches the other side, where 
it will turn back into an electron. This phenomenon is 
called Klein tunneling. 

An explanation for this phenomenon is that increas- 
ing barrier height leads to an increased degree of mode- 
matching between the wavefunctions of the holes within 
the barrier and the electrons outside of it. When the 
modes are perfectly matched (in the case of an infinitely 
tall barrier), we have perfect transmission through the 
barrier. In the case of graphene, the chirality discussed 
earlier leads to a varying transmission probability de- 



Some experimental results have been interpreted as ev- 
idence for Klein tunneling. Klein tunneling has been ob- 
served through electrostatic barriers, which were created 
by gate voltages (Stander et al. 2009). Similar effects 



have also been observed in narrow graphene resonant het- 



erostructures (Young and Kim 2009). 



D. Graphene vs traditional materials 

Here we summarize some of the interesting properties 
of graphene by comparing them with more traditional 
materials such as 2D semiconductors. 

1. Traditional semiconductors have a finite band gap 
while graphene has a nominal gap of zero. Nor- 
mally, the study of electron and hole motion 
through a semiconductor must be done with dif- 
ferently doped materials. However, in graphene 
the nature of a charge carrier changes at the Dirac 
point from an electron to a hole or vice- versa. On 
a related note, the Fermi level in graphene is al- 
ways within the conduction or valence band while 
in traditional semiconductors the Fermi level often 
falls within the band gap when pinned by impurity 
states. 

2. Dispersion in graphene is chiral. This is related to 
some very distinctive material behaviors like Klein 
tunneling. 

3. Graphene has a linear dispersion relation while 
semiconductors tend to have quadratic dispersion. 
Many of the impressive physical and electronic 
properties of graphene can be considered to be con- 
sequences of this fact. 

4. Graphene is much thinner than a traditional 2D 
electron gas (2DEG). A traditional 2DEG in a 
quantum well or heterostructure tends to have an 
effective thickness around 5-50 nm. This is due to 
the constraints on construction and the fact that 
the confined electron wavefunctions have an evanes- 
cent tail that stretches into the barriers. Graphene 
on the other hand is only a single layer of carbon 
atoms, generally regarded to have a thickness of 
about 3 A (twice the carbon-carbon bond length). 
Electrons conducting through graphene are con- 
strained in the z-axis to a much greater extent than 
those that conduct through a traditional 2DEG. 



5. Graphene has been found to have a finite minimum 
conductivity, even in the case of vanishing charge 
carriers dNovoselov et a/.[ |2005| |Tan et a/.[|2007[ ). 
This is an issue for the construction of field-effect 
transistors (FETs), as we will see in more detail 
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in section |VI| since it contributes to relatively low 
on/off ratios for graphene-based transistors. 

Readers seeking further reading about distinctly elec- 
tronic properties of graphene should refer to Sec. |VI[ The 
next section will introduce and discuss the vibrational 
properties of graphene. 



III. VIBRATIONAL PROPERTIES 

While the electronic properties have attracted the 
lion's share of the interest in graphene, the vibrational 
properties are of great importance too. They are re- 
sponsible for several fascinating properties such as record 
thermal conductivities. Since graphene is composed of a 
light atom, where the in-plane bonding is very strong, 
graphene exhibits a very high sound velocity. This large 
sound velocity is responsible for the very high thermal 
conductivity of graphene that is useful in many appli- 
cations. Moreover, vibrational properties are instrumen- 
tal in understanding other graphene attributes, including 
optical properties via phonon-photon scattering (e.g. in 
Raman scattering) and electronic properties via electron- 
phonon scattering. 



A. Phonon dispersion 

Most of the vibrational properties of graphene can be 
understood with the help of the phonon dispersion re- 
lation. Interestingly, the phonon dispersion has some 
similarity with the electronic band structure discussed 
in the previous section, which stems from the identical 
honeycomb structure (the out-of-plane modes are shown 
in Fig. [3|. In order to obtain the phonon dispersion it is 
necessary to consider the vibrational modes of the crys- 
tal in thermal equilibrium. This is done by considering 
the displacement of each atom from its equilibrium po- 
sition, written Un for the atom labeled n. Each atom 
is effectively coupled to its neighbors by some torsional 
and longitudinal force constants, which only depend on 
the relative positions of the atoms. This allows one to 
write the Newtonian coupled equation of motion in fre- 
quency space as: 



(4) 



where in graphene the sum over m is typically over 
the second or fourth nearest neighbors, with the corre- 
sponding coefficients ^rn,n, also known as the dynamical 
matrix. 

In graphene, and similarly to the electronic structure, 
the two sublattices A and B have to be considered ex- 
plicitly to solve for the eigenspectrum of the dynamical 



matrix. However, the atoms can vibrate in all three di- 
mensions, hence the dynamical matrix has to be written 
in terms of both the sublattices A and B as well as the 3 
spatial dimensions. This leads to a dynamical matrix in 
reciprocal space ^rn,n{k)^ which is given by a 6 x 6 ma- 
trix when assuming u^'^ ~ e*^'^^'^. Here one applies 
Bloch's theorem, where R^'^ is the equilibrium position 
of atom n in the sublattice A and B, respectively. Two of 
the eigenvalues correspond to the out-of-plane vibrations, 
ZA (acoustic) and ZO (optical), and the remaining 4 cor- 
respond to the in-plane vibrations: TA (transverse acous- 
tic), TO (transverse optical), LA (longitudinal acoustic) 
and LO (longitudinal optical). 
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FIG. 3 First Brillouin zone and out-of-plane phonon 
modes. The phonon dispersion relation of the ZA and ZO 
modes as a function of the in-plane reciprocal vector k. The 
vertical axis is the phonon frequency, while the horizontal axes 
are momentum space on the graphene lattice. The dispersion 
relation is obtain ed using the second-nearest-neighbor model 
for graphene ( Falkovsky 2007 ) . The gray surface corresponds 



to the ZO (optical) mode, whereas the pink surface shows the 
ZA (acoustic) mode. Also shown are the corresponding K, F, 
and M points of the Brillouin zone. 

The ZA and ZO modes are often assumed to be decou- 



pled from the in-plane modes ( Falkovsky[ |2007), which 



leads to a simple dispersion relation very similar to the 
electronic band structure as shown in Fig. [3) It is inter- 
esting to note that the dispersion is quadratic at the F 
point, which is unusual for acoustic modes. In contrast, 
a simple graphene phonon model based on atomic poten- 
tials containing only three parameters leads to a linear 



dispersion for the ZA mode (Adamyan and Zavalniuk 



2011 ). However, the experimental data seems to be more 
consistent with a quadratic dispersion (see Fig.[5| ( Dres- 
selhaus and Eklund[ |2QQQ[ |Popov[ |2QQ4| ). At the K and 
K' points we recover a cone structure similar to the Dirac 
cones in the electronic structure. However, the phonon 
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density of states does not vanish at these points because 
of the presence of the in-plane modes. 
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FIG. 4 First Brillouin zone and in-plane phonon 
modes. The phonon dispersion relation of the TO and LO 
modes in gray and the TA and LA modes in pink as a function 
of the in-plane reciprocal vector k. The longitudinal modes 
are on top of the transverse modes. The vertical axis is the 
phonon frequency, while the horizontal axes are the momen- 
tum space on the graphene lattice. 

The in-plane modes are constituted by two acoustic 
modes and two optical modes. These modes can be 
obtained from the reduced in-plane dynamical matrix, 
which can be described by a 4 x 4 matrix, assuming no 
coupling to the out-of-plane modes. Using the param- 



eters given by reference (Falkovsky, 2007), the in-plane 
modes are shown in Fig. [4] 




FIG. 5 Experimental phonon dispersion relation in 
graphite. All the phonon modes are shown, including the 
in-plane and out-of-plane modes. The data (full symbols) is 
shown along the special symmetry points F, K, and M along 
with r esults from ab initio calc ulations (lines). Figure adapted 
from (iWirtz and Rubiol 12004). 



The in-plane modes show the expected linear disper- 
sion at the F point. The transverse modes closely fol- 



low the longitudinal modes but with a slightly lower fre- 
quency, for both the acoustic and optical modes. For 
comparison with experiments, it is more instructive to 
show the dispersion relation following surface cuts along 
the lines F to M, M to K and K to F in the Brillouin 
zone. An extensive collection of experimental data for 
graphite is shown in Fig. [5] along with ab initio calcula- 
tions. The data was obtained using several techniques, 
including neutron scattering, electron energy loss spec- 
troscopy. X-ray scattering, infrared absorption and dou- 



ble resonant Raman scattering experiments (Wirtz and 



Rubio , 2004 ) . There is no comparably extensive data on 
graphene yet, but it is expected to be very similar to 
graphite, since the coupling between planes is very weak. 
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FIG. 6 Theoretical phonon dispersion of graphene and 
Raman spectrum of large scale graphene. The labels 
(G, D, 2D, etc) identify the peaks from the Raman spectrum 



with the corresponding phonon energies (Rao et al. 2011 



Venezuela et al. \ |2011| ) . Only the in-plane modes are shown 



since they are the only ones which are Raman active. Several 
Raman data sets of graphene ^"^C are shown in blue taken 
from different spots of the same sample at a laser wavelength 
of 514nm. The red line corresponds to the average over these 
data sets. The phonon dispersion is calculated using the mass 
of ^^C. 



Most of the data on graphene stems from Raman scat- 
tering, which allows for the determination of the phonon 
spectrum close to special symmetry points. This is illus- 
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trated in Fig. |6j where the theoretical phonon dispersion 
with parameters from (Falkovsky, 2007) is shown with 
the labels corresponding to the various Raman peaks. 
The Raman peaks determine the phonon energies for 
some values of the momentum. The Raman data was ob- 
tained from chemical vapor deposition (CVD) ^^C grown 
graphene. More conventional ^^C graphene is very simi- 
lar except for a rescaling of the Raman peaks due to the 
change of mass. Other aspects of the Raman spectrum 
are discussed in Sec. [Vl 

An important consequence of the phonon dispersion re- 
lation in graphene is the very high value of the in-plane 
sound velocity, close to Cph 20 km/s (Adamyan and 



Zavalniuk 2011), which leads to very high thermal con- 



ductivities. 



B. Thermal conductivity 

From the kinetic theory of gases, the thermal con- 
ductivity due to phonons is given hy ^ Cph Cv{T) A, 
where (7y(T) is the specific heat per unit volume and A 
is the phonon mean free path. This implies that since 
Cph is very large in graphene, one can expect a large 
thermal conductivity. Indeed, experiments at near room 
temperature obtain n 3080-5150 W/mK and a phonon 
mean free path of A 775 nm for a set of graphene 
fiakes dBalandin efZ] |2008| [Ghosh et a/l|2008D . 



These results indicate that graphene is a good can- 
didate for applications to electronic devices, since a 
high thermal conductivity facilitates the diffusion of heat 
to the contacts and allows for more compact circuits. 
Phonons also play an important role in electronic trans- 
port via electron-phonon scattering, which is discussed 



in Sec. VI The mechanical properties of graphene are 



discussed in Sec. |VIII[ whereas the synthesis of graphene 
is treated in the next section. 



the cost and difficulty involved in the method. 

Table |l] shows a summary of some of the most impor- 
tant synthesis methods. The typical number of graphene 
layers produced as well as currently achievable dimen- 
sions are given. For comparison the mobilities listed are 
for graphene transferred to Si/Si02 wafers, since the elec- 
tron mobility of graphene is heavily substrate dependent, 
as discussed in Sec. [Vll 



Method 


Layers 


Size 


Mobility (cm^V-^s-^) 


Exfoliation 


1 to 10+ 


1 mm" 


15000^ 


Thermal SiC 


1 to 4 


50 /im^ 


2000^ 


Ni-CVD 


1 to 4 


1 cm^ 


3700^ 


Cu-CVD 


1 


65 cm^ 


16000^ 



TABLE I Comparison of graphene synthesis methods. 

Shows typical number of layers produced, size of graphene 



laye rs (larg est dimension) and mo bilit y on Si/Si02. a : ( Geim 
|2009i), b:(|Novoselov a/.| |2005|, c:(|Emtsev et al 



d:( |Kim a/.[ |2009| ), e:( |Bae et a/.l|2Qlb| ), f:( |Li et a/.| (2010^ . 
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A. Mechanical exfoliation 

Developed by Geim and Novoselov, the exfoliation pro- 
cess uses HOPG (highly oriented pyrolitic graphite) as 
a precursor. The HOPG was subjected to an oxygen 
plasma etching to create 5 jiia deep mesas and these 
mesas were then pressed into a layer of photoresist. The 
photoresist was baked and the HOPG was cleaved from 
the resist. Scotch tape was used to repeatedly peel flakes 
of graphite from the mesas. These thin flakes were then 
released in acetone and captured on the surface of a 



Si/Si02 wafer (Novoselov et a/., 2004) 



IV. SYNTHESIS 

Since graphene was isolated in 2004 by Geim and 
Novoselov using the now famous Scotch tape method, 
there have been many processes developed to produce 
few-to-single layer graphene. One of the primary con- 
cerns in graphene synthesis is producing samples with 
high carrier mobility and low density of defects. To date 
there is no method that can match mechanical exfoliation 
for producing high-quality, high- mobility graphene flakes. 
However, mechanical exfoliation is a time consuming pro- 
cess limited to small scale production. There is great in- 
terest in producing large scale graphene suitable for ap- 
plications in flexible transparent electronics, transistors, 
etc. Some concerns in producing large scale graphene are 
the quality and consistency between samples as well as 




FIG. 7 Monolayer graphene produced by mechanical 
exfoliation. Large sample with length of 1 mm on Si/Si02. 

These few-layer graphene (FLG) flakes were identified 
using the contrast difference in an optical microscope and 
single layers using an SEM. Using this technique Geim 
and Novoselov were able to generate few- and single- 
layer graphene flakes with dimensions of up to 10 /im 
(Novoselov et al., 2004). The few-layer graphene flakes 



were found to have ballistic transport at room temper- 
ature and mobilities as high as 15000 cm^V~^s~^ on 
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Si/Si02 wafers (Novoselov et al, 2005). The scotch tape 
method can generate flakes with sides of up to 1 mm in 
length ( Geim 2009 ) , of excehent quahty and weh suited 
for fundamental research. However, the process is lim- 
ited to small sizes and cannot be scaled for industrial 
production. 



B. Thermal decomposition of SiC 

The thermal decomposition of silicon carbide is a tech- 
nique that consists of heating SiC in ultra-high vacuum 
(UHV) to temperatures between 1000°C and 1500°C. 
This causes Si to sublimate from the material and leave 
behind a carbon rich surface. Low-energy electron mi- 
croscopy (LEEM) studies indicate that this carbon layer 
is graphitic in nature, which suggests that the technique 



could be used to form graphene (Hass et al. 2008|). 



Berger and De Heer produced few-layer graphene by 
thermal decomposition of SiC. The Si face of a 6H-SiC 
single crystal was first prepared by oxidation or H2 etch- 
ing in order to improve surface quality. The sample was 
then heated by electron bombardment in UHV to 1000° C 
to remove the oxide layer. Once the oxide was removed 
the samples were heated to 1250-1450° C, resulting in the 
formation of thin graphitic layers. Typically between 1 
and 3 layers were formed depending on the decomposition 
temperature. Using this method, dev ices were produce d 
with mobilities of 1100 cm V-^s-^ (Berger et a/.l|2004|). 



This technique is capable of generating wafer-scale 
graphene layers and is potentially of interest to the semi- 
conductor industry. Several issues still remain, notably 
controlling the number of layers produced, repeatability 
of large area growths and interface effects with the SiC 



substrate (Choi et al. 2010). 



Emtsev et al. found that by heating SiC in Ar at 900 
mbar as opposed to UHV they were able to reduce surface 
roughness and produce much larger continuous graphene 
layers, up to 50 jiia in length. The graphene on SiC 
was characterized using atomic force microscopy (AFM) 
and LEEM, as shown in Fig. [Sj Th ey measured electro n 



mobilities of up to 2000 cm^V-^s"^ ( [Emtsev et a/.[|2009 ). 

Juang et al. synthesized millimeter size few- to single- 
layer graphene sheets using SiC substrate coated in a thin 
Ni film. 200 nm of Ni was evaporated onto the surface of 
the SiC and the sample was heated to 750° C in vacuum. 
Graphene was found to segregate to the surface of the Ni 
on cooling. This gave a continuous graphene layer over 



2009) 



the entire nickel surface (Juang et al. 

Unarunotai et al. have developed a technique to trans- 
fer graphene synthesized on SiC onto arbitrary insulating 
substrates. Graphene was first produced using a typical 
thermal decomposition of SiC technique. A bilayer film 
of gold/polyimide was deposited onto the SiC wafer and 
then peeled off. The gold/polyimide film was then trans- 
ferred onto a Si/Si02 substrate and the gold/polyimide 





FIG. 8 Graphene produced by thermal decomposition 

of SiC. (a) AFM image of graphene growth on SiC annealed 
at UHV. (b) LEEM image of UHV grown graphene film, (c) 
AFM image of graphene annealed in Ar at 900 mbar. (d) 
LEEM image of graphene on Ar ann ealed SiC substrate s how- 
ing terraces up to 50 /xm in length ( Emtsev et al. 2009|) . 



layers were removed using oxygen plasma reactive ion 
etching. This yielded single-layer graphene fiakes with 



mm areas ( Unarunotai et al. , 2009 ) 



C. Chemical vapor deposition 

In contrast to the thermal decomposition of SiC, where 
carbon is already present in the substrate, in chemical 
vapor deposition (CVD), carbon is supplied in gas form 
and a metal is used as both catalyst and substrate to 
grow the graphene layer. 



1. Growth on nickel 

Yu et al. grew few-layer graphene sheets on polycrys- 
talline Ni foils. The foils were first annealed in hydrogen 
and then exposed to a CH4-Ar-H2 environment at atmo- 
spheric pressure for 20 mins at a temperature of 1000° C. 
The foils were then cooled at different rates between 
20°C/s and 0.1°C/s. The thickness of the graphene layers 
was found to be dependent on the cooling rate, with few- 
layer graphene (typically 3-4 layers) being produced with 
a cooling rate of 10°C/s. Faster cooling rates result in 
thicker graphite layers, whereas slower cooling prevents 
carbon from segregating to the surface of the Ni foil (Yu 



et 



20086 



To transfer the graphene layers to an insulating sub- 
strate, the Ni foil with graphene was first coated in sil- 
icone rubber and covered with a glass slide then the Ni 
was etched in HNO3. 
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2. Growth on copper 

Li et al. used a similar process to produce large scale 
monolayer graphene on copper foils. 25 fim thick copper 
foils were first heated to 1000° C in a flow of 2 seem (stan- 
dard cubic centimeters per minute) hydrogen at low pres- 
sure and then exposed to methane flow of 35 seem and 
pressure of 500 mTorr. Raman spectroscopy and SEM 
imaging confirm the graphene to be primarily monolayer 
independent of growth time. This indicates that the pro- 
cess is surface mediated and self limiting. They fabri- 
cated dual gated FETs using graphene and extracted a 



carrier mobility of 4050 cm^y-^s-^ (iLi et a/.M2009). 



Recently a roll-to-roll process was demonstrated to 
produce graphene layers with a diagonal of up to 30 
inches as well as transfer them to transparent flexible 
substrates (Bae et al 2010). Graphene was grown by 



CVD on copper, and a polymer support layer was ad- 
hered to the graphene-copper. The copper was then re- 
moved by chemical etching and the graphene film trans- 
ferred to a polyethylene terephthalate (PET) substrate. 
These films demonstrate excellent sheet resistances of 125 
r^/n for a single layer. Using a repeated transfer process, 
doped 4-layer graphene sheets were produced with sheet 
resistances as low as 30 Q/B and optical transmittance 
greater than 90%. These 4-layer graphene sheets are su- 
perior to commercially available indium tin oxide (ITO) 
currently used in flat panel displays and touch screens in 
terms of sheet resistance (~100 Q/B for ITO) and optical 
transmittance (-90% for ITO). 




FIG. 9 Multiple CVD graphene sheets transferred to 

PET. A roll-to-roll process was used to produce graphene 



sheets with up to 30 inch diagonal (Bae et al. 2010) 



Li et al. have shown the dependence on the size of 
graphene domains synthesized by CVD with tempera- 
ture, methane flow and methane pressure. Performing 
the growth at 1035° C with methane flow of 7 seem and 
pressure 160 mTorr led to the largest graphene domains 
with average areas of 142 /im^. A two-step process was 



used to first grow large graphene flakes and then by mod- 
ifying the growth conditions to fill in the gaps in the 
graphene sheet. Using this technique they were able to 
produce sa mples with carrie r mobility of up to 16000 
cm^V~^s~^ ('Li et al. 2010b). In general, the graphene 



layer is slightly strained on the copper foil due to the 



high temperature growth (Yu et al. 2011). 




FIG. 10 Controlling domain size in CVD graphene. Ef- 
fect of temperature, methane flow and methane partial pres- 
sure on the size of graphene dom ains in CVD growth, scale 
bars are 10 /xm (|Li et al. 2010b). 



Recently Lee et al. have demonstrated a technique 
to produce uniform bilayer graphene by chemical vapor 
deposition on copper using a similar process but with 
modified growth conditions. They determined optimal 
bilayer growth conditions to be: 15 minutes at 1000° C 
with methane flow of 70 seem and pressure of 500 mTorr. 
The bilayer nature of the graphene was confirmed by Ra- 
man spectroscopy, AFM, and transmission electron mi- 
croscopy (TEM). Electrical transport measurements in a 
dual gated device indicate that a band gap is opened in 



CVD bilayer graphene (Lee et al. 2010). 
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FIG. 11 Bilayer CVD growth on copper, (a) 2 x 2 inch 
bilayer graphene on Si/Si02. (b) Raman spectrum with 514 
nm laser source of 1 and 2 layers of graphene produced by 
exfohation and CVD (iLee et a/.l 12010). 



F. Sodium-ethanol pyrolysis 

Graphene was produced by heating sodium and 
ethanol at a 1:1 molar ratio in a sealed vessel. The prod- 
uct of this reaction is then pyrolized to produce a material 
consisting of fused graphene sheets, which can then be re- 
leased by sonication. This yielded graphene sheets with 
dimensions of up to 10 jiia. The individual layer, crys- 
talline and graphitic nature of the samples was confirmed 
by TEM, selected area electron diffraction (SAED) and 



Raman spectroscopy (Choucair et a/., 2009). 



G. Other methods 



D. Molecular beam deposition 

Zhan et al. succeeded in layer-by-layer growth of 
graphene using a molecular beam deposition technique. 
Starting with an ethylene gas source, gas was broken 
down at 1200° C using a thermal cracker and deposited 
on a nickel substrate. Large area, high quality graphene 
layers were produced at 800° G. This technique is capable 
of forming one layer on top of another, allowing for syn- 
thesis of one to several layers of graphene. The number 
of graphene layers produced was found to be independent 
of cooling rate, indicating that carbon was not absorbed 
into the bulk of the Ni as in GVD growth on nickel. 
Results were confirmed using Raman spectroscopy and 



TEM (Zhan et a/., 2011) 
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FIG. 12 Molecular beam deposition produced 
graphene. (a) Diagram of thermal crack er setup, (b ) TEM 
image of graphene film, scale bar 100 nm (Zhan et al. 2011). 



E. Unzipping carbon nanotubes 

Multi- walled carbon nanotubes were cut longitudinally 
by first suspending them in sulphuric acid and then treat- 
ing them with KMn04. This produced oxidized graphene 
nanoribbons which were subsequently reduced chemi- 
cally. The resulting graphene nanoribbons were found 
to be conducting, but electronically inferior to large scale 
graphene sheets due to the presence of oxygen defect sites 



(Kosynkin et al 2009). 



There are several other ways to produce graphene such 
as electron beam-irradiation of PMM A nanofibres ( Duan 



et al. 



2008), arc discharge of graphite ( Subrahmanyam 



et al. , 2009), thermal fusion of PAHs (Wang et al. , 2008a), 



and conversion of nanodiamond (Subrahmanyam et al. 



2008). 



H. Graphene oxide 

Another approach to the production of graphene is son- 
ication and reduction of graphene oxide (GO). The polar 
O and OH groups formed during the oxidation process 
render graphite oxide hydrophilic, and it can be chemi- 
cally exfoliated in several solvents, including water ( Zhu| 
al, [2010,). The graphite oxide solution can then be 
sonicated in order to form GO nanoplatelets. The oxy- 
gen groups can then be removed in a reduction process 
involving one of several reducing agents. This method 
was used by Stankovich et al. using a hydrazine reduc- 
ing agent, but the reduction process was found to be 
incomplete, leaving some oxygen remaining (Stankovich] 



et 



2007 ). 

Graphene oxide (GO) is produced as a precursor to 
graphene synthesis. GO is useful because its individual 
layers are hydrophilic, in contrast to graphite. GO is sus- 



pended in water by sonication (McAllister et al, 2007 



^Paredes et al. 2008), then deposited onto surfaces by 
spin coating or filtration to make single or double layer 
graphene oxide. Graphene films are then made by re- 
ducing the graphene oxide either thermally or chemically 



(Marcano et al., 2010). The exact structure of graphene 
oxide is still a matter of debate, although there is consid- 
erable agreement as to the general types and proportion 



of oxygen bonds present in the graphene lattice ( He et al. 



1996). 



1. Wet chemical synthesis 

The chemical methods to produce GO were all devel- 
oped before 1960. The most recent and most commonly 
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employed is the Hummers procedure (Hummers and Offe- 2. Plasma functionalization 



man 



1958 ). This process treats graphite in an anhydrous 



mixture of sulfuric acid, sodium nitrate, and potassium 
permanganate for several hours, followed by the addi- 
tion of water. The resulting material is graphite oxide 
hydrate, which contains approximately 23% water. Sub- 
sequent nuclear magnetic resonance and X-ray diffraction 
studies of the structure of GO have led to fairly detailed 
models based on a combination of hydroxide, carbonyl, 
carboxyl and epoxide groups covalently bonded to the 
graphene lattice. Fig. 13 shows a predicted structure 



for GO produced using the Hummers method (He et al. 



1998). 



graphite oxide 




FIG. 13 Structure of a monolayer of graphite oxide (He 
et a/.| 119981). 



Following the realization of the potential importance 
of graphene as a replacement for semiconductor materi- 
als and indium tin oxide (ITO), as discussed in Sec. X.B 



alternative methods for graphene production have been 
explored. Other approaches have been sought to produce 
the same hydrophilicity in graphite without the time and 
material requirements of the Hummers method. Very re- 
cently, glow discharge treatment has been proven to in- 
troduce oxygen species into the lattices of all forms of 
graphitic materials (e.g. buckyballs, CNTs, graphene, 
carbon nanofibers, and graphite) ( [Vandsburger et aT. 
2009). The resulting graphene/graphite oxides have a 



structure very similar to Hummers GO, and can be ther- 
mally treated to selectively reduce epoxides. Unlike the 
Hummers method, plasma functionalization requires no 
strong acids, can proceed at room temperature and can 
be completed very quickly, often in a matter of seconds 
or minutes. 

Aside from its potential to replace the Hummers 
method, plasma treatment in itself is interesting for al- 
tering the electrical conductivity of graphene or thin 
graphite. This allows for bandgap engineering as well 
as phenomena like photoluminescence (PL). 
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FIG. 14 Summary of the Hummers method and ther- 
mal reduction. Bulk graphite is oxidized then separated 
in water. Then it is thermally reduced to make single- layer 
graphene ([McAllister et a/. | [2007). 



3. RF plasma 

Radio frequency (RF) plasma refers to a processing 
technique whereby a capacitive plasma is ignited in an 
isolated volume. RF treatment is used almost exclusively 
for surface treatment of graphene, because ion bombard- 
ment is significantly reduced in RF treatments, as op- 
posed to DC discharges. In RF plasma, electrodes need 
not be in contact with the plasma gas, current is sup- 
plied with an alternation frequency of 13.57 MHz, and 
power ranges from 10 W to 50 W. RF treatment has 
been shown to selectively affect the outermost surface of 
graphene (Hazra et al. 2011). Other work using only 



oxygen for RF treatment allows for layer-by-layer etch- 



ing of a graphite surface, producing islands of GO (You 



et al 1993). 



Understandably, the degree of oxidation strongly af- 
fects the in-plane electrical and thermal conductivity of 
graphene oxide. Increased introduction of oxygen groups 
into the graphene lattice interrupts the sp'^ hybridiza- 
tion of electron orbitals. Epoxide groups can be reduced 
by thermal treatment, or reaction with potassium iodide 
(KI), leading to a similar structure, where only hydroxyl 
groups are present. This leads to improved electrical con- 
ductivity and unaffected hydrophilicity. In each case the 
procedure requires tight temperature control and long 
reaction times of several hours. The basic process, in- 



cluding thermal treatment, is shown in Fig. 14 



4. Photoluminescence 

Single- and dual-layer graphene do not exhibit photo- 
luminescence, due mainly to the negligible bandgap of 
native graphene. GO does have a photoluminescent re- 
sponse, but the typical oxidation methods, sonication of 
bulk graphite oxide, are inappropriate for use in photo- 
luminescence applications. For that reason, RF plasma 
oxidation has been the subject of recent work at produc- 
ing photoluminescent single layers of GO. The procedure 
for producing GO thin films from single layer graphene 
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was reported in (Gokus et al, 2009). Rather than oxi- 



dizing bulk graphite to produce single layers of graphene, 
single- or few-layer graphene is oxidized after isolation. 
Typically, graphene is prepared by micro-cleavage using 
the scotch tape or other methods, and electrical contacts 
or other additions are put in place. Following installa- 
tion, RF plasma treatment in Ar-02 mixes is applied in 
one to six second intervals. The plasma power reported 
was 10 W, the pressure 0.04 mbar, and the gas composi- 
tion ratio was 2:1 Ar to O2. 



An important finding reported by (Childres et al. 



2011) is the temporal evolution of the Raman spectrum 
of graphene with increasing plasma treatment time. The 
spectra are reproduced in Fig. [TsJ The most noticeable 
change in the spectra is the gradual reduction in inten- 
sity of the 2D and 2D' peaks, which are indicative of 
sp^ hybridization. This indicates the disruption of the 
graphene lattice by introduction of oxygen groups and 
demonstrates oxidation. Further changes in the region of 
interest involve the development of a G peak that arises 
from the increased presence of "disordered" carbon. Sim- 
ilar findings were reported in other work using pulsed RF 
plasma treatment rather than continuous treatments. 
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FIG. 15 Raman spectra of graphene w ith increasing 
num ber of 1.5 sec plasma treatments ( Childres et al] 

2oTT|. 



In a sample with a heterogeneous surface, it was shown 
that only regions of single-layer graphene oxide displayed 
PL, while untreated graphene or multi-layer graphene did 
not. The first image shown in Fig. 16 is an image of 
the PL produced by laser fiuorescence. The bright re- 
gions correspond to low intensity regions in an elastic 
scattering image (b), revealing that they are single-layer 
graphene, labeled IL in the chart, (c) shows both a PL 
curve and a contrast curve taken along the white dotted 



line in (a). The middle-contrast regions in the blue curve 
correspond exactly to high PL. 
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FIG. 16 Photoluminescence in oxygen plasma treated 

GO. (a) Dark regions are pristine graphene, while bright re- 
gions are single layer graphene oxide, (b) Few-layer graphene 
is bright while single-layer graphene oxide is dark, (c) Shaded 
zone shows the correlation between contrast and PL-intensity 
taken along the white dotted line in (a) (Gokus et al. 2009). 



Photoluminescence occurs after plasma treatment as 
a result of the introduction of defects in the graphene 
lattice. Such defects disrupt the electrical properties 
of pristine graphene and introduce a bandgap that is 
absent from native graphene monolayers. A bandgap is 
desirable for more than PL, and other work has reported 
plasma oxidation of single and few layer graphene for 
these purposes (Nourbakhsh et al, 2010). Specifically, 



a bandgap would allow for logic and optoelectronic 



applications as we shall see respectively in Sec. [IX] and 
Sec.Kl 



Once graphene has been produced, it is important to 
identify it and to characterize its structure, which is the 
topic of the next section. 
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V. CHARACTERIZATION 



A great many techniques are being used to characterize 
graphene. We discuss here some of the most important 
ones with a particular emphasis on the identification of 
graphene. 



A. Raman spectroscopy 

Raman spectroscopy is an important characterization 
tool used to probe the phonon spectrum of graphene 
as discussed in section Sec. |III[ Raman spectroscopy 
of graphene can be used to determine the number of 
graphene layers and stacking order as well as density of 
defects and impurities. The three most prominent peaks 
in the Raman spectrum of graphene and other graphitic 
materials are the G band at ~1580 cm~^, the 2D band at 
~2680 cm~^ and the disorder-induced D band at ~1350 
cm~^. 

The G band results from in-plane vibration of sp^ car- 
bon atoms and is the most prominent feature of most 
graphitic materials. This resonance corresponds to the 
in-plane optical phonons at the T point. The 2D band 
arises as a result of a two phonon resonance process, in- 
volving phonons near the K point, and is very promi- 



nent in graphene as compared to bulk graphite ( Ni et al 



2008a|). 

The D band is induced by defects in the graphene lat- 
tice (corresponding to the in-plane optical phonons near 
the K point), and is not seen in highly ordered graphene 
layers. The intensity ratio of the G and D band can be 
used to characterize the number of defects in a graphene 



sample (Pimenta et al, 2007). 



The line shape of the 2D peak, as well as its intensity 
relative to the G peak, can be used to characterize the 
number of layers of graphene present as illustrated in 
Fig. [TtI Single-layer graphene is characterized by a very 
sharp, symmetric, Lorentzian 2D peak with an intensity 
greater than twice the G peak. As the number of layers 
increases the 2D peak becomes broader, less symmetric 



and decreases in intensity (Wang et al. 2008b). 



B. Optical microscopy 

Monolayer graphene becomes visible on Si02 using an 
optical microscope. The contrast depends on the thick- 



ness of Si02, the wavelength of light used ( Blake et al. 
[2007^ and the angle of illumination ( ,Yu and Hilke[|2009D ". 



This feature of graphene is useful for the quick identifi- 
cation of few- to single-layer graphene sheets, and is very 
important for mechanical exfoliation. Fig. [18] shows the 
optical contrast of one, two and three layers of exfoliated 
graphene under different wavelengths of illumination and 
different thicknesses of Si02. 
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FIG. 17 Layer dependence of graphene Raman spec- 
trum. Raman spectra of N = 1-4 layers of graphene on 
Si/Si 02 and of bulk graphite. Figure adapted from ^Yu[ 
20T0l. 
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FIG. 18 Optical microscope images of graphene. Mul- 
tilayer graphene sheet on Si/Si02 showing optical c ontrast at 
different wavelengths and thicknesses (Blake et al. 2QQ7J . 



C. Electron microscopy 



Transmission electron microscopy has been used to im- 
age single-layer graphene suspended on a microfabricated 
scaffold. It was found that single-layer graphene dis- 
played long range crystalline order despite the lack of 



a supporting substrate ( [Meyer et al. , 2007). Suspended 
graphene was found to have considerable surface rough- 
ness with out-of-plane deformations of up to 1 nm. 

Aberration-corrected annular dark-field scanning 
transmission electron microscopy (ADF-STEM) was 
used in order to image CVD grown graphene suspended 
on a TEM grid ( [Huang et Z| [2011[ ). They found that 



along grain boundaries the hexagonal lattice structure 
breaks down and the grains are "stitched together" with 



pent agon- heptagon pairs as seen in Fig. 20 
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FIG. 19 Atomic scale TEM image of suspended 
graphene. Few- to single-layer graphene she et show ing long 
range crystalline order, scale bar 1 nm (Meyer et aL\ (2007). 
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FIG. 20 ADF-STEM imaging of graphene suspended 

on TEM grid, (a) SEM image of graphene transferred to 
TEM grid, scale bar 5 /xm. (b) Atomic scale ADF-STEM im- 
age showing the hexagonal lattice in the interior of a graphene 
grain, (c) ADF-STEM image showing intersection of two 
grains with a relative rotation of 27°. (d) Same image with 
pentagons, heptagons and deformed hexagons formed along 
grai n boundary highlig hted (b, c and d have scale bars of 5 
A) |Huang et ar||2011 ). 



D. Measuring the electronic band structure 

A wide variety of experimental techniques exist for 
measuring the band structure of materials. Due to its 
particular characteristics, graphene places severe limita- 
tions on the techniques available. Most band structure 
measurement techniques are highly sensitive to the bulk 
of a material rather than the surface. Since graphene 
is so thin, we need techniques that are very sensitive to 
surface layers. 

Angle-resolved photoemission spectroscopy (ARPES) 



is the most popular technique for measuring the band 
structure of graphene. Photons of sufficient energy (20- 
100 eV) are shot at the surface of the material being 
probed. Each photon is energetic enough that if it inter- 
acts with an electron, it has a significant chance of trans- 
ferring enough energy to launch the electron out of the 
material completely. The electron must be given enough 
energy to overcome the work function of the material. 

The electron, once free of the material, will have a 
chance of hitting the ARPES detector. The detector 
is oriented so that it can measure one specific angle of 
electron emission. Note that there are two degrees of 
freedom in angle, typically called (j) and 0. Using these 
two together, one can specify any direction. The detec- 
tor is also able to accurately measure the energy E of 
the outgoing electron. This means that ARPES will si- 
multaneously measure the three variables 0, ^, and E. 
The three components (x,y,z) of the scattered electron's 
momentum prior to being struck by the photon can be 
found using the measured quantities. In this way the 
experimenter can map out the correspondence between 
energy and momentum within the material with high res- 
olution. 




FIG. 21 Band structure of graphene on top of SiC. 

Vertical axis is the electron's energy, and horizontal axis is 
its momentum. Note the key locations in momentum space 
with reference to Fig. [2] and also that F is zone center, corre- 
sponding to zero momentum. The black line is a theoretical 
prediction based on the tight binding approximation. The 
fainter bands are believed to be due to interactions between 
the substrate and the graphene. Image adapted from (Bost- 



wick et al] 2007) 



ARPES is capable of scanning to within about 5 A 
of the surface when using electrons of 20-100 eV. This 
means that most of the signal will be from the first 
few atomic layers of the surface in question. This prop- 
erty makes ARPES particularly well-suited to measuring 
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FIG. 22 Substrate-induced band gap in single layer 
graphene on top of SiC. (a) Real space and momentum 
space structure of graphene. (b) Band structure of graphene 
taken along vertical black line near the K point in panel (a). 
The black lines are dispersion relations estimated from energy 
distribution curves. Figure from ( ,Zhou et al 2007). 



the band structure of incredibly thin materials such as 
graphene. On the other hand, this introduces some ex- 
perimental difficulties because it means that the sample 
surface must be kept under ultra-high vacuum (UHV). 
Creating and measuring graphene without leaving UHV 
is a significant experimental challenge. Another com- 
monly used method is to anneal the graphene by running 
a significant electrical current through it. The annealing 
process does a good job of cleaning the graphene so that 
it can be measured by techniques like ARPES even after 
it has been exposed to atmosphere. 

ARPES measurements have been made on graphene in 
a wide variety of circumstances and with many different 
goals. For the purposes of this review, only a few studies 
will be mentioned out of this vast field. Fig. [21] shows the 
experimental band structure of graphene grown on top of 



SiC ( [Bostwick et al. , 2007). The intent of this study was 
to better understand the dynamics of charge carriers in 



graphene. Fig. 22 is from a different group that was also 
probing the behavior of graphene grown on top of SiC 



(Zhou et a/., 2007). This second experiment observed a 
notable band gap in their single-layer graphene samples. 
Additionally, they noticed that this gap shrank as the 
number of layers of graphene was increased from one to 
four. It is believed that the existence of the observed 
band gap is due to interactions with the substrate that 
cause the symmetry of graphene's 7r-bonds to be broken 



(Zhou et a/., 2008). 



There are a number of variations of ARPES that 
differ only in the wavelength of the probing photons. 
Angle-resolved ultraviolet photoemission spectroscopy 
(ARUPS) has also been used to study graphene's band 
structure ( [Gierz et al. 2008). The primary reason why 
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TABLE II Comparison between the electronic prop- 
erties of graphene and common bulk semiconductors. 

Energy band gap (Eg), electron effective mass (m* /me), elec- 
tron mobility (/ie) and electron saturation velocity (lygat) of 
graphene is compared to those of conventional semiconductors 
and AlGaN/GaN 2DEG ( |Giannazzo et a/.||2011| . 



this technique is employed is convenience. In ARUPS, 
a laboratory-based ultraviolet wave source can be used 
to produce the probing photons. This is a less expensive 
and simpler setup than ARPES, which typically uses X- 
rays produced from a synchrotron. It is also worth not- 
ing that information about the electronic structure of 
graphene can be inferred from the results of other mate- 



rial techniques such as optical spectroscopy (Mak et al. 
20101). 



In the previous sections we have discussed various ways 
to obtain and identify graphene. We now turn our atten- 
tion to its physical properties, starting with electronic 
transport measurements. 



VI. ELECTRONIC TRANSPORT AND FIELD EFFECT 

Owing to its unique band structure (see Sec. [ill ), 
graphene exhibits novel transport effects such as ambipo- 
lar field effect and minimum conductivity which are ab- 



sent in most conventional materials (Wu et al. 2010). 
This unusual electronic behavior leads to exceptional 
transport properties in comparison to common semicon- 
ductors. This can be seen on Table [ll] which compares 
two of the main electronic properties (carrier mobility 
and saturated velocity) of graphene with those of com- 
mon bulk semiconductors and 2DEGs. In what follows, 
we will first describe the experimental methods that are 
commonly used to measure the ambipolar field effect. We 
will then discuss the transport properties that can be ex- 
tracted from this experimental data. The effect of dif- 
ferent scattering mechanisms on the carrier mobility and 
minimum conductivity will then be discussed in detail. 
Finally, other electrical properties relevant to transistor 
technology will be reviewed. 
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A. Measurement of the ambipolar field effect 

Transport properties are typically measured with a 
graphene device similar to those shown in Fig. |23j To 
fabricate these devices, graphene (exfoliated, CVD, etc) 
is often deposited on an oxidized silicon wafer (Si02/Si). 
Later on we discuss other substrates that are sometimes 



used (Dean a/. 2010b Ponomarenko a/. 2009). Un- 



less otherwise specified, all measurements reported here 
were made using exfoliated graphene flakes. Electrical 
contacts, usually made of gold, are then deflned using a 
lithographic process or a stencil mask to avoid photore- 
sist contamination. Electrodes are generally patterned 
in a 4-lead (Fig. [23^ ) or Hall bar (Fig. [23|3) conflgura- 
tion. Lastly, the device can be cleaned by annealing at 
ultrahigh vacuum or in H2/Ar gas, or by applying a large 
current density (~ 10^ A/cm^) through it to remove ad- 
sorbed contamination ( |Moser et al. 2007). 





FIG. 23 Schematic representation of common elec- 
tronic d evices, (a) 4-lead (Dorgan et al. 2010J and (b) 
Hall bar ( [Novoselov et al.\ |2004| ). 



With this graphene device in hand, one can tune the 
charge carrier density between holes and electrons by ap- 
plying a gate voltage {Vg) between the (doped) silicon 
substrate and the graphene flake. The gate voltage in- 
duces a surface charge density n = e^eVg/te where e^e 
is the the permittivity of Si02, e is the electron charge 
and t is the thickness of the Si02 layer. This charge den- 
sity change shifts accordingly the Fermi level position 
{Ef) in the band structure (see the insets of Fig. [24^). 
At the Dirac point, n should theoretically vanish (Gein? 



and Novoselov [ 2007), but as will be explained further 
on, thermally generated carriers {nth) and electrostatic 
spatial inhomogeneity (n*) limit the minimum charge 
density (i Dorgan et al\ |2010 ). Fig. 24 3, which shows 
the calculated carrier density as a function of gate volt- 
age, clearly illustrates the fact that charge density is well 
controlled by the gate away from the Dirac point. The 
linear relation between Vg and n was verifled experimen- 
tally ( Novoselov et al. 2004[ ) in that region by measuring 
the Hall coefficient Rh = 1/ne as a function of Vg (see 
Sec. VII. A.2). Typically, charge density can be tuned 
from 10'' to 10^^ cm by applying a gate voltage that 



moves Ef 10 to 400 meV away from the Dirac point (Gi- 



annazzo et al. 2011 ). 



(a) 
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1 1 1 



(b) 2 




FIG. 24 Ambipolar electric field effect in graphene 

The insets of (a) show the changes in th e position of the Fermi 
level ^/ as a function of gate voltage (Geim and Novoselov 



2007|. (b) Cal culated charge density vs. gate voltage at 300 
K and 500 K (Dorgan et al. |2010"). Solid lines include con- 
tribution from n*, nth and Ug. Dashed line shows only the 
contribution from the gating (rig). 



In 2004, the ambipolar fleld effect corresponding to 
the change in resistivity p (or conductivity a = 1/p) 
that occurs when the charge density is modified by the 



gate voltage was observed and analyzed ( Novoselov et al. 



2004). Experimentally, p is measured using a standard 
4-probe technique ( [Van Der Pauw| |1958[ ) and is given 
by p = {W/ L){V2^/ Iia) where W and L are respectively 
the width and the length of graphene, is the volt- 
age across electrodes 2 and 3 (see Fig. 23) and /14 is 
the current between contact 1 and 4. Note that because 
of the uncertainty on the aspect ratio L/W^ the error 
on the absolute magnitude of p is usually around 10% 
( Chen et a/. , 2008a|b ). As Fig. 24^ shows, resistivity 
rapidly increases as we remove charge carriers, reaching 
its maximum value at the Dirac point. From this curve 
one can extract the carrier mobility p = 1/enp and the 
minimum conductivity (Jmin' Other definitions for mo- 
bility are sometimes used such as the field effect mobility 
fJ^FE = {l/C)da/dVg (where C is the gate capacitance) 
( Dean et al.\ 201Qbl) an d the Hall mobility pHaii = Rh / P 
( Novoselov et a/.[|2004 ). Note that in practice, the carrier 
mobility is only meaningful away from the Dirac point, 
where n is accurately tuned by the gate voltage. 



B. Transport and scattering mechanisms 



In contrast with the ideal, theoretical graphene, exper- 



and ripples (Katsnelson and Geim 2008) and is affected 



by phonons (Bolotin et al. 2008a). These perturbations 



imental graphene contains defects ( Chen et al. , 2009 ) and 
impurities (Chen et al., 2008a Zhang et al. , |2009a ), in- 
teracts with the substrate (Chen et al. , 2008b), has edges 



alter the electronic properties of a perfect graphene sheet 
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first by introducing spatial inhomogeneities in the carrier 
density and, second, by acting as scattering sources which 



reduce the electron mean free path (Giannazzo et al 



2011 ). The former effect dominates when the Fermi level 



is close to the Dirac point and alters the minimum con- 
ductivity of graphene whereas the latter effect prevails 
away from the Dirac point and affects the carrier mo- 
bility. The impact of these perturbations has been sub- 
jected to intensive and ongoing investigation, on both the 
theoretical and experimental side, in order to determine 
the mechanisms that limit the mobility and the mini- 
mum conductivity. From a theoretical point of view, two 
transport regimes are often considered depending on the 
mean free path length / and the graphene length L. When 
/ > L, transport is said to be ballistic since carriers can 
travel at Fermi velocity {uf) from one electrode to the 
other without scattering. In this regime, transport is de- 



scribed by the Landauer formalism ( Peres , 2009 ) and the 
conductivity is expressed as: 



1. Phonon scattering 

Phonons can be considered an intrinsic scattering 
source since they limit the mobility at finite temperature 
even when there is no extrinsic scatterer. As explained in 
section Sec. the dispersion relation of graphene com- 
prises six branches. Longitudinal acoustic (LA) phonons 
are known to have a higher electron-phonon scattering 
cross-section than those in the other branches (| Hwang 
and Das Sarma, 2008). The scattering of electrons by LA 



phonons can be considered quasi-elastic since the phonon 
energies huq are negligible in comparison with Ep, the 
Fermi energy of electrons. 

In order to determine the effect of electron-phonon 
scattering on resistivity, one must consider two distinct 
transport regimes separated by a characteristic temper- 
ature Tbg called the Block- Griineissen temperature, de- 



fined as ( Hwang and Das Sarma , 2008 ) : 



CThall 



(5) 



where the sum is over all available transport modes 
of transmission probability T^. For ballistic transport 
mediated by evanescent modes, this theory predicts that 
at the Dirac point the minimum conductivity is: 



4e^ 



4.92 X 10"^ Q-^ 



(6) 



On the other hand, when / < L, carriers undergo elas- 
tic and inelastic collisions and transport enters the diffu- 
sive regime. This regime prevails when the carrier den- 
sity n is much larger than the impurity density n^. In 
that case, transport is often described by the semiclassi- 
cal Boltzmann transport theory (Das Sarma et al. 2011) 



and at very low temperature carrier mobility can be ex- 
pressed in terms of the total relaxation time r as: 



(8) 



where /c^ is the Boltzmann constant, Vph is the sound 
velocity and kp is the Fermi wave vector with reference 
to the K point in the BZ. 



(9) 



where n is the electron density in the conduction 



band (Pisana et al. 

cm~^ we get Tbg 



n 



10 



2007|). If one measures n in units of 
54v/n K. 



Consider first T ^ TbGi the equipartition {EP) limit. 
In this case the Bose-Einstein distribution function for 
the phonons is N{uOq) ^ kBT/hUq^ which leads to a lin- 
ear T-dependence of the scattering rate and hence the 
resistivity p ^ T. In the BG or degenerate regime, on 
the other ha nd, where T <^ Tbg^ one obtain s at very low 
p ^ (Hwang and Das Sarma 2008) as shown in 
Fig.|25l 



h 



(7) 2. Coulomb scattering 



This equation describes the diffusive motion of carriers 
scattering independently off various impurities. The re- 
laxation time depends on the scattering mechanism dom- 
inating the carrier transport or a combination thereof. 
The scattering mechanisms mostly discussed in the lit- 
erature include Coulomb scattering by charged impuri- 
ties (long range scattering), short-range scattering (de- 
fects, adsorbates) and electron-phonon scattering. In the 
following, we provide a brief theoretical introduction of 
these scattering processes and relevant transport mea- 
surements. 



Coulomb scattering stems from long-ranged variations 
in the electrostatic potential caused by the presence of 
charged impurities close to the graphene sheet. These im- 
purities are often thought of as trapped ions in the under- 
lying substrate and screened by the conduction electrons 
of graphene. Assuming random distribution of charged 
impurities with density Ui and employing a semiclassical 



approach, it was predicted ( [Adam et al] |2007j that the 
charged-impurity scattering is proportional to ^ . With 
Eq. ([t]), the conductivity at high carrier density (n ^ n^) 
is given by: 
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T(K) 



FIG. 25 Electric resistivity of graphene at ultrahigh 
carrier densities. Resistivity over wide range of T, showing 
the cross-over from the low T ^ (BG ) regime to the high 
T- (EP) one ( |Efetov and Kim[|201Q| . 
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(10) 



where C is a dimensionless parameter related to the 
scattering strength. Considering the random phase 
approximation and the dielectric screening from the 
Si02 substrate, it was predicted that C ~ 20. 

Chen et al. investigated experimentally the effect of 
charged impurities on the carrier mobility and conduc- 
tivity by doping a graphene flake with a controlled potas- 



sium flux in UHV ( Chen et al 2008a ) . Fig. 26 i shows the 



conductivity as a function of gate voltage for a pristine 
sample and three different doping concentrations. It can 
be clearly seen that the gate voltage of minimum conduc- 
tivity becomes more negative with increasing doping. As 
it was previously shown (Schedin et a/.[|2007a[ ), this is be- 
cause K atoms dope graphene with electrons (n-doping), 
which in effect moves the Fermi level up with respect to 
the Dirac point. From Fig. [26^ , one can also see that 
(jiVg) becomes more linear and mobility decreases as the 
doping concentration increases which is in good agree- 



ment with Eq. (10). The dashed line in Fig. 26 3 shows 



that mobility scales linearly with l/n^ when transport 
is limited by charged-impurity scattering. In these mea- 
surements C ~ 20 in Eq. (10) was obtained (Chen et al. 



2008a Tan et al ^ 2007). 



FIG. 26 Effect of charge impurities and defects on 
transport properties of graphene. (a) The conductivity 
(cr) vs. gate voltage for a pristine sample (black curve) and 
three different potassium doping concentrations at 20 K in 
UHV. Lines represent empirical fits (Chen et a/., 2008a|). (b) 
Inverse of mobility (l//i) as a function of ion dosage for differ- 
ent samples and irradiated ions (Ne+ and He+). The dashed 
line corresponds to t he behavior for the same concentration 
of potassium doping (Chen et al 2009). 



3. Short-range scattering 

Finally, short-range defects such as vacancies and 
cracks in graphene flakes are predicted to produce 
midgap states in graphene (Stauber et al, 2007). Va- 



cancies can be modeled as a deep circular potential well 
of radius R and this strong disorder gives rise to a con- 
ductivity which is roughly linear in n: 



(11) 



where rid is the defect density. This equation mim- 



ics the one for charged impurities (Eq. (10)), with a 
slightly logarithmic dependence of the conductivity on 
the charge carrier density. Defect scattering was exper- 
imentally studied by irradiating a graphene flake with 



500 eV He and Ne ions in UHV (Chen et al, 2009). The 



resulting conductivity was also demonstrated to be ap- 
proximately linear with charge density, with mobility in- 
versely proportional to the ion (or defect) dose Ud. As 
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shown in Fig. [26p , the mobiUty decrease was found to 
be 4 times larger than the same concentration of charged 
impurities. From the hnear fits of Fig. 26 with Eq. (11), 
the impurity radius was found to be ~ 2.9 A which is 
a reasonable value for a single-carbon vacancy. 



C. Mobility 

1. Graphene on Si02 

In a graphene electronic device, all of the scattering 
mechanisms mentioned above come into play. From a 
technological point of view, determining the exact nature 
of the scattering that limits the mobility is essential 
in order to develop high-speed electronic devices. To 
do so, one must also take into account the effect of 
the underlying substrate on the electronic transport. 
Chen et al performed a general study of scattering 
mechanisms in graphene on Si02 by measuring the 



temperature dependence of the ambipolar effect ( |Chen 



et 



2QQ8b[ ). The measurements were fitted using 
three terms: p{n,T) = po{n) + Pla{T) + pspp{n,T). 
Each term was associated with a certain scattering 



mechanism. As shown in Fig. 27 i, the first two terms 
were determined with a linear fit at low T. The 
y-intercepts correspond to po{n) and were found to 
scale linearly with . According to Eq. (10) and 



Eq. (11), this behavior can be associated with charged 
impurities and defects. Pla{T) was determined by 
the slope of the linear fit in Fig. |27| and appears to 
be independent of the charge density. These results 
agree very well with predictions for longitudinal acoustic 
(LA) phonon scattering in the regime where T ^ Tbg- 
Finally, the third term pspp^ was found to have a 
strong temperature dependence. Such behavior was 
explained by considering the surface polar phonons 
(SPP) of the Si02 substrate which produce an electrical 
field that couples to electrons in graphene. Theoreti- 
cal expressions for SPP-limited resistivity showed good 



agreement with the data of Fig. 27 i at high temperature. 



Fig. [27p presents the temperature-dependent mobility 
and the theoretical limits of the three scattering mecha- 
nisms (LA phonons, SPP and impurities/defects). One 
can see that at room temperature, SPP and impuri- 
ties/defects are by far the two dominant scattering mech- 
anisms for graphene on Si02. At low and room tem- 
perature, mobility is mainly limited by the former type 
of scattering. Typically, mobility of graphene on Si02 
ranges from 10000 to 15000 cm^V~^s~^. It was suggested 
that the scattering impurities are trapped charges in the 
underlying Si02 substrate, but this remains under debate 
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FIG. 27 Temperature dependence of resistivity and 
mobility in graphene on Si02. (a) Temperature- 
dependent resistivity of graphene on Si02 for different gate 
voltages, samples and runs. Dashed lines are linear fits 
at low temperature, (b) Temperature-dependent mobility in 
graphene on Si02 for two samples at n = 1 x 10^^ cm~^. Mo- 
bility limits due to LA phonons (dark red solid line), substrate 
surface phonons (green dashed line), and impurities/defects 
(red and blue dashed lines) are presented. Matthiessen's rule 
was used to obta in the net mobility for each sample (red and 
blue solid lines) (|Chen et a/.[ |2008b|). 



2. Suspended graphene 

Removing the substrate or using one which is free of 
trapped charge are two possible ways to improve the car- 
rier mobility. The former approach was used by (Du 
et a/.[|2QQ8 ), in which the authors fabricated a suspended 
graphene device by chemically etching the underlying 
Si02 (see the inset of Fig. 28 i). From this device, they 



obtained mobilities as high as 200,000 cm V s for 



charge density below 5x10 cm . Fig. 28 i shows the 



( Ponomarenko et al. , 2009 ) . 



charge density dependence of mobility for suspended and 
non-suspended (on Si02) graphene devices at 100 K. The 
arrows indicate the minimum charge density that is ac- 
curately controlled by the gate voltage. From this fig- 
ure, one can see that mobility in suspended graphene 
approaches the ballistic model prediction. Furthermore, 
the mobility deterioration caused by adsorbed impurities 
on the surface of suspended graphene has been demon- 
strated (Bolotin et a/., 2008a|b ). As Fig. 28 3 shows, mo- 
bility and mean free path / clearly increase after current- 
induced cleaning and high-temperature annealing of the 
device. At large carrier density (n > 0.5 x 10^^ cm~^) and 
temperatures above 50 K, resistivity of clean suspended 
graphene was found to be linear with temperature, sug- 
gesting LA phonon scattering. At low temperature (~ 5 
K), mobility as high as 170,000 cm^V~^s~^ was obtained 
and the mean free path reached the device dimension. In 
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FIG. 28 Conductance and mobility in suspended 
graphene as a function of charge density, (a) Mobility 
vs. charge density for suspended (red line) and non-suspended 
(black line) graphene at T = 100 K. The blue line represents 
the ballistic model prediction. Inset: schematic representa- 
tion of the suspended graphene device ( Du et al.^ 2008). (b) 
Conductance vs. gate voltage (at T = 40 K) for a suspended 
graphene device before (blue line) and after (red line) anneal- 
ing and current-induced cleaning. The red dotted line was 
calculated using a ballistic mo del. Inset: AFM im age of the 
suspended graphene device ( |Bolotin et al. 2008a). 



graphene/BN devices are troublesome to fabricate and 
are thus not ideal for industrial applications. Currently, 
large area graphene synthesized by CVD or thermal seg- 
regation is the most promising material for technological 
applications since it can be produced on a wafer scale 
dKim et al\ |2009| |Li et al\ |2009D . Depending on the 



technique, large scale graphene (on Si02) typically show s 
mobilities lower than 5000 cm^V~^s~^ (Avouris 2010). 



The scattering mechanisms responsible for this reduced 
mobility are still under investigation. Preliminary stud- 



ies suggested that grain boundaries (Yazyev and Louie 



2010), doping defects (Wei et al. 



gate voltage control (Cao et al. 



2009) and ineffective 



2010) might limit the 



mobility. In addition to mobility, a common way to char- 
acterize the quality of large scale graphene is to measure 
its sheet resistance Rg (resistivity at = 0). This prop- 
erty varies with the synthesizing technique used and can 
be as low as 125 ^t/U ( |Bae et aT\ 2010). It was also shown 
that sheet resistance increases with decreasing tempera- 



ture ( [Park et al. , 2010). This observation contrasts with 
the metallic behavior {da/dT < 0) of pristine graphene 
flakes but agrees with the insulating behavior of irradi- 



ated graphene (with defects) (Chen et a/., 2009). 



Substrate Production 
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TABLE III Mobility range (/i) and minimum con- 
ductivity (amin) of graphene produced by different 
techniques and deposited on different substrates, a: 


(Novoselov et al. 


2004), b: 


(Dean et al. 2010b), c: (Bolotin 



f: i\U et a/.|f2010b| ) 



these conditions, conductivity is well described by the 
ballistic model as shown in Fig. [28b. 



D. Minimum conductivity 



3. Other substrates 

Although suspended graphene shows impressive trans- 
port properties, this geometry imposes evident con- 
straints on the device architecture. To overcome this 
problem, boron nitride (BN) was proposed as a sub- 
strate ( [Dean et aJ] |201Qb[ ). Compared to Si02, BN 
has an atomically smooth surface, is relatively free of 
charged impurities, has a lattice constant similar to 
that of graphene and high surface phonon frequencies. 
All these advantages result in a mobility about three 
times higher than that of graphene on Si02. However, 



The presence of disorder in graphene (ripples, defects, 
impurities, etc.) produces fluctuations in its electro- 
static potential. We mentioned above that these per- 
turbations are significant at the Dirac point where the 
screening of the potential fluctuations is weak due to 
the low charge density. These fluctuations in the charge 
density can be thought of as electron-hole puddles and 
have been observed experimentally using scanning probe 



methods (Martin et al.^ 2008 Zhang et al.^ 2009a) on 



graphene/Si02 samples. Experimentally, these disor- 
dered graphene samples have a minimum conductivity 
about TT times larger than that predicted for ballistic 
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transport given in Eq. (|6|. This discrepancy between 
theory and experiments was observed in 2004 ( [Novoselov 
et al , 2004 ) and was known as "the mystery of the miss- 
ing pi." The effect of doping on the minimum conductiv- 
ity was also investigated in a study of charged-impurity 
scattering ( Chen et al. , 2008a ) . It was found that (Jmin 



decreases on initial doping and reaches a minimum near 
Ae^/h only for non-zero charged impurities. This sug- 
gests that charged impurities trapped in the Si02 sub- 
strate, located at the graphene/Si02 interface or on the 
graphene surface are responsible for the minimum con- 
ductivity obtained experimentally. 

To verify this conclusion, the minimum conductivity 
of clean suspended graphene was measured as a function 
of temperature (Du et a/., 2008|). As Fig. 29 i shows. 



minimum conductivity decreases with temperature and 
approaches the ballistic prediction down to a factor of 
1.7. However, the linear temperature dependence was 
not expected. Another group also investigated the near- 
ballistic regime by measuring the minimum conductivity 
on samples with different aspect ratios (W/L) and surface 
areas at low temperature (~1.5 K) ( |Miao et al] [2007 ). 
For large- area samples (~3 /im^), the minimum conduc- 
tivity is mostly independent of the aspect ratio (inset of 
Fig. 293). On the other hand, small-area devices (<0.2 



/im^) yield results that are clearly dependent on the as- 
pect ratio (Fig. [29|3). Devices with wide electrodes and 
short channels (large W/L) approach the theoretical min- 
imum conductivity in the ballistic regime. Table |TIT| com- 
pares the carrier mobility and minimum conductivity of 
graphene produced by different techniques and on differ- 
ent substrates. 



E. Other transport properties 

Concerning electronic applications, graphene has at- 
tracted considerable attention due to its high mobility. 
However, other transport properties must be taken into 
account in order to develop graphene-based technologies 
successfully. Among those properties, saturation velocity 
is particularly relevant for field-effect transistor (FET) 
applications. In modern FETs, short channel lengths re- 
sult in high electrical fields ( [Schwierz 2010) of around 
~100 kV/cm. In such conditions, the carriers acquire 
enough kinetic energy to excite the optical phonon modes 



of graphene (Avouris 2010). As a consequence, carrier 



velocity saturates, reducing the relevance of mobility to 
device performance. The optical phonons of graphene 
have higher energy (~160 meV) than those of common 
semiconductors such as Si (55 meV) ( [Dorgan et al. , 2010). 
Consequently, the intrinsic saturation velocity is higher 
in graphene than in conventional semiconductors. 

Numerical values of saturation velocity are presented 



in Table [Hj and Fig. [30] shows the experimental and the- 
oretical dependence of saturation velocity on charge den- 
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FIG. 29 Minimum conductivity of graphene. (a) Mini- 
mum conductivity vs. temperature for a suspended graphene 
device. The upper dashed line represents the 4e^/h limit for 
graphene with charged impurities. The lower dashed line cor- 
responds to the theoretical value in the ballistic regime (Du 
\et al.\ |2008). (b) Minimum conductivity vs. as pect ratio for 



small and large (inset) graphene/Si02 devices (Miao et al. 
2007t . 



sity for graphene/Si02. A simple phonon emission model 
(Meric et al. |2011| ) predicts that saturation velocity is 
proportional to In Fig. 30, the upper and lower 

theoretical curves correspond to saturation velocities lim- 
ited by optical phonons of graphene and Si02, respec- 
tively. The experimental results suggest that both kinds 
of phonons play a role in limiting Ugat but that substrate 



phonons of Si02 are dominant for this device ( ^Dorgan 
et a/., 2010). Some studies also point out that current 



never reaches a complete saturation in some graphene de- 
vices. This incomplete saturation might be due to a com- 
petition between disorder and optical phonon scattering 



(Barreiro et a/., 2009), or the formation of a "pinch-off" 



region (Meric et a/., 2011). Mechanisms limiting the sat- 
uration velocity are under active investigation. 

Graphene has excellent transport properties, but to 
take advantage of them, carriers must be injected and col- 
lected through metal contacts. These electrical contacts 
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FIG. 30 Electron saturation velocity of graphene on 

Si02. Electron saturation velocity vs. charge density for an 
electric field of 2 V//im at T = 80 K and 300 K. The upper 
and lower dashed curves correspond to theoretical saturation 
velocities limi ted by optical phon ons of graphene and Si02, 
respectively ( (Dorgan et al.\ 2010). 



VII. MAGNETORESISTANCE AND QUANTUM HALL 
EFFECT 

The first demonstration, in the early 1980s, of the ex- 
istence of a so-called quantum Hall effect (QHE) of the 
two-dimensional electron gas (2DEG) at the interfaces 
of semiconductor heterostructures sparked great inter- 
est in the experimental study of the properties of low- 
dimensional systems (Klitzing et a/., 1980). Thus it is 



perhaps no surprise that when the one atom thick, two- 
dimensional lattice of graphene could finally be isolated 
in 2004, physicists immediately tried to characterize its 
magnetotransport properties, which are still an impor- 
tant field of research today The purpose of this sec- 
tion is to recapitulate the present state of experimen- 
tal knowledge about magnetotransport measurements in 
graphene. Firstly, we will describe the general experi- 
mental procedures that are necessary to carry such mea- 
surements. Secondly, we will briefly review the funda- 
mental physics relevant to localization phenomena and to 
the QHE. Finally, we will show how these techniques pro- 
vide information on the electronic properties of graphene 
through weak (anti-) localization, the integer QHE and 
the fractional quantum Hall effect (FQHE). 



A. Experimental procedures 



produce energy barriers that limit the charge transfer at 
the graphene/metal junction f Avouris , 2010). This limi- 
tation results in a contact resistance which can be mea- 
sured and should be reduced to create high-performance 
graphene devices. One way to determine the contact re- 
sistance is to fabricate variable channel length devices 
and to extrapolate the resistance to zero channel length. 
It was shown that contact resistance is temperature and 
gate voltage dependent, and that it varies from about 
100 ft /im to a few kQ /im (jXia et al 2011). Finally, 
it is worth pointing out that graphene can sustain cur- 
rent densities greater than 10^ A/cm^, which is 100 times 



higher than those supported by copper (Moser et al 
[2007). Graphene can thus be used as interconnects in 
integrated circuits. 



Since Geim's and Novoselov's seminal work (Novoselov 



et al , 2004), graphene's unique electronic properties have 



attracted massive interest and created an explosion of 
scientific activity. Much ink has been spilled about the 
ambipolar field effect, the ultra-high mobility of graphene 
and the limiting scattering mechanisms, as well as the 
minimum conductivity. In this section, we reported some 
of the main experimental studies on these ever growing 
subjects. In the next section, we turn our attention to 
the related but more specific topics of magnetotransport 
and quantum Hall effect in graphene. 



We first describe the basic experimental procedures 
required to perform low temperature magnetotransport 
measurements in graphene. 



1. Setup 

The setup consists of a superconducting electromagnet 
(typically Nb based type-II superconductor, e.g. NbTi or 
NbSna) immersed in a liquid helium cryostat. It is rep- 



resented in Fig. 31 Thus the sample can be refrigerated 
to temperatures of 4 K and below. Lower temperatures 
are achieved by reducing the vapor pressure of helium by 
pumping. Furthermore, in liquid helium the fields gener- 
ated in the electromagnet can reach high maximal values 
typically ranging between 6 and 14 T. 

The sample is placed at the bottom of the cryostat in 
such a way that the field is perpendicular to the plane 
of the graphene. The device is usually conceived for 4- 



terminal measurements, as illustrated in Fig. 32 A small 
current is driven through the end terminals and voltage 
measurements are carried out between different combina- 
tions of the lateral terminals. This avoids the measure- 
ment of the contact resistance at the current injection 
points. The longitudinal resistivity pxx is the ratio of 
the longitudinal voltage to the current normalized by the 
aspect ratio. The Hall resistivity pxy and the Hall Re- 
sistance Rxy are identical and are defined as the ratio 
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graphene sheet is placed on a Si02 substrate about 330 
nm thick, which in turn hes on a doped Si wafer. A gate 
voltage is applied on the Si to inject electron carriers in 
the graphene or withdraw them. In the latter case, holes 
carry the current. Thus graphene can be studied above 
and below the neutral point corresponding to the Fermi 
energy and for which in principle n = 0. 

The calibration is done using the conventional Hall ef- 
fect. In a (low) applied magnetic field, the Hall resis- 
tance is related to the current by Rxy = RrBI^ where 
Rh = is the Hall constant and q is the charge of the 
carriers. A measurement of for different gate volt- 
ages Vg allows one to establish a linear map between the 
carrier density n and Vg for both holes and electrons. An 



example of calibration is shown on Fig. 33 



FIG. 31 Setup for magnetotransport measurements. 

The cryostat is cooled to T = 4 K in two steps using liq- 
uid nitrogen and liquid helium. A solenoidal superconducting 
magnet creates a magnetic field perpendicular to the plane of 
the graphene. The helium bath can be pumped on to yield 
even lower temperatures. 



of the transverse voltage to the current. The Hall con- 
ductivity cFxy is related to the resistivity by the inverse 
tensor relation cFj.y = 2 - 




FIG. 32 Magnetotransport measurements. Top: 
Graphene is etched into a 4-terminal Hall bar. Gold electrodes 
are used to perform electrical measurements. The width of 
the central wire is 0.2 /im. Bottom: Graphene is placed on 
a thin Si02 layer (285-300 nm), which is itself on a doped 
Si wafer. A gate vo ltage Vg allows the car rier density to be 
tuned. Figure from (Novoselov et al. 20051. 



2. Carrier density tuning 

In graphene, the carrier density n can be tuned us- 
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FIG. 33 Measurement of the conventional Hall effect 

for different gate voltages Vg gives Vg (x.n. The Hall constant 
Rh is related to the carrier density by Rh — Figure 
from (Novoselov et al. 2005). 



3. Observation of the quantum Hall effect 

The observation of the QHE is conditional on many 
factors. 

1. High quality samples must be used in order to max- 
imize the momentum scattering time Tq of the elec- 
trons. 

2. High magnetic fields B are required in order for 
the cyclotron period of the carriers to be much 
shorter than the scattering time Tq. Equivalently, 
this means that the cyclotron frequency ujc of an 
electron must be much higher than the broadening 
of the carrier energy levels (the Landau levels, see 
Sec.|VII.B[): cjcTq > 1. 



ing the electric field effect, as discussed in Sec. VI A 



3. Low temperatures are needed in order for the ther- 
mal energy ksT to be much less than the spacing 
hujc of the Landau levels: ksT << hujc- 
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In graphene, the integer QHE can be observed at liq- 
uid helium temperature for fields ranging from about 1 
T and higher. For high magnetic fields the effect can be 



closed trajectory (Castro Neto et al, 2009). This is key 



observed for much higher temperatures ( [Novoselov et al. 
2005 ). As an extreme example it has been shown that for 



to understanding the weak localization measurements in 
graphene discussed below. 

Since in graphene there are two inequivalent Dirac val- 
leys K and K' with the same dispersion, there is an equa- 



a gigantic field of 5 = 45 T, the QHE remains detectable tion similar to Eq. ( 12 ) for K'. This leads to an extra val 



at room temperature (Novoselov et al. , 2007). The obser- 



vation of the FQHE requires much more extreme condi- 
tions. The FQHE can be observed in graphene for fields 
above 2 T for temperatures lower than 1 K and it can sub- 
sist up to T = 20 K at B = 12 T, provided that the sam- 
ple is ultraclean and that the graphene is suspended to 
suppress scattering. This yields very low charge inhomo- 
geneity of order n = 10^ cm~^ (Du et al , i2009). In such 



small devices the contacts of the 4-terminal measure- 
ment prevent the appearance of the QHE and one must 



perform a two-terminal measurement instead ( |Du et al. 



2009). More recently, FQHE measurements in graphene 



ley degeneracy for the state of the electron and to chiral 
electrons: the state in valley K has a different he- 
licity than a state ^j^/ in valley K'. We call the doublet 
(K, K') the valley pseudospin. We conclude that the en- 
ergy levels should be 4 times degenerate: two times for 
the pseudospin and two times for the actual spin of the 
electron. 

In a graphene bilayer, two sheets of graphene are 
stacked onto each other (in the natural Bernal stacking 
of graphite) (Castro Neto et al. 2009). In such a system, 
the carriers can hop between the layers and the dispersion 
relation is not linear. Instead, the valence and conduction 
bands consist of two parabolic bands of the same curva- 



have been performed on graphene lying on a hexagonal 
boron-nitride substrate at fields of the order of B = 35 

T and temperatures of the order of T = 0.3 K (|Dean 2009). The theory of Eq. (12) can then be extended to 



ture touching at the neutral point (Castro Neto et al. 



et al., 2010a). Such fields cannot be reached with the 
simple setup illustrated in Fig. [31] Instead, these mea- 



surements are performed in High Magnetic Field Labora- 
tories such as the National High Magnetic Field Labora- 
tory at Florida State University or Grenoble High Mag- 
netic Field Laboratory in France. 



B. Theoretical background 

1. Chiral electrons and pseudospin 



show that the carriers should be massive Dirac fermions 



with a Berry phase of 2ti ( ,Novoselov et al. 2006[ ) 



2. Landau levels 



Eq. (12) can be solved in the presence of a magnetic 



field and the result predicts a sequence of energy levels 
called the Landau levels (LL) ([Apalkov and Chakraborty 



20061 [Castro Neto et a/.||2009| |Ezawa||2008||Gusynin and| 
Sharapovl |2005[ ): 



Imagine a gas of electrons in the x-y plane confined by 
some potential well in the z-direction. We assume that 
this well is deep and narrow enough, and hence that the 
energy levels of the well are separated enough, so that 
the electrons are forbidden to access its excited states 
through any excitations. Then we say that the electron 
motion is two-dimensional. 

At low carrier density n, the dispersion of graphene is 
linear around a Dirac point (call it K) so that we can- 
not define an effective mass such that E = as in 

2m* 

semiconductor samples. Instead it can be shown that in 
the continuum limit (i.e. taking the position of the elec- 
tron on the lattice to be a continuous variable), the mo- 
tion of the electron around a single Dirac point should 
be described by the Dirac- Weyl equations for massless 
fermions discussed in Sec. [ll| usi ng the Hamiltonian of 
Eq. (|3| dCastro Neto et a/.[]2Q09l ): 



(12) 



where ^ k is a spinor. It can be shown that the state 
^ K acquires an extra phase (a Berry phase) of tt on a 



En = sgn{N)^2envl \N\ TV = 0, ±1, ±2... (13) 

The energy is measured with respect to the Dirac point 
(i.e. the energy at the Dirac point is zero). Each such 
level should also have the 4-fold degeneracy discussed 
above. Note that in a semiconductor 2DEG, the LLs 
would be equally separated by an amount TiuJc where ujc = 
is the cyclotron frequency and rUc the cyclotron mass 



(Datta, 1997), which from Eq. (13) is clearly not the case 



for graphene. We sketch a figure of the splitting of the 



bands in LLs for graphene in Fig. 34 



For bilayer graphene, we find a different expression: 



En = sgn{N)nuJc^\N\{\N\^l) N 



0,±1,±2... 

(14) 



It can be shown that the extra layer degree of freedom 
should give rise to an extra degeneracy for the level = 
and only for that level. Thus, the A" = level is 8 times 
degenerate while the others keep the 4-fold degeneracy of 
graphene. The splitting of the bands in LLs for bilayer 
graphene is illustrated schematically in Fig. l34| 
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FIG. 34 Schematic representation of the formation of 
Landau levels in graphene (top) and bilayer graphene (bot- 
tom) . In both cases the levels are not equally spaced and the 
neutral point has an associated LL. In graphene each level is 
4 times degenerate because of the spin and valley pseudospin. 
In bilayer graphene the lowest LL A/" = is 8 times degenerate 
because of the extra layer degree of freedom. 



C. Measurements of weak localization and antilocalization 
in graphene 

1. Quantum interferences 

Quantum interferences have long been known to affect 
transport measurements at low temperatures. In partic- 
ular, two-dimensional electron systems in the presence 
of a (low) magnetic field show variations of conductance 
with respect to the value in the absence of a field. This 
effect is known as weak (anti-) localization. 

Consider an electron hitting an impurity and going 
back in the direction from which it came. In a classical 
picture, this could happen in a number of ways. In partic- 
ular, the electron could go around the impurity clockwise 



or counter-clockwise (see Fig. 36). However, since elec- 
trons are quantum mechanical waves, we must add each 
of the possible paths and make them interfere together to 
get the net probability of the electron to have backscat- 
tered. Electrons in conventional 2DEGs, in the absence 
of spin-orbit interaction and scattering by magnetic im- 
purities, gain the same phase on both trajectories and in- 
terfere constructively, leading to an increased probability 
of backscattering compared to the value expected from 
the Drude model. We say that the electrons exhibit weak 
localization (WL). If a magnetic field is applied, an addi- 
tional (Aharonov-Bohm) phase is added between the two 
paths, destructive interference occurs and the conduc- 
tance increases (positive magnetoconductance) (Tikho- 



nenko et 



2009t 



In graphene, however, it can be shown that because 
of the additional Berry phase of tt of the wave func- 
tion, the two paths end up in opposite phase in the ab- 
sence of spin-orbit interaction and scattering within or 
between Dirac valleys. Therefore application of a mag- 
netic field can only restore the constructive interference 
and decrease the conductance (negative magnetoconduc- 



tance) ( [Morozov et a/.[ |2QQ6| [Tikhonenko et a/.[|2QQ9 ). 

We say that the electrons exhibit weak anti-localization 
(WAL). Note that conventional electrons which couple 
strongly through spin-orbit interactions can also exhibit 
WAL. However, this effect is not expected in graphene 



because of the small mass of the carbon atom (Tikho- 



nenko et al 2009). 



2. Measurements 

The early quantum interference measurements on 
graphene showed no WAL and a strongly suppressed WL. 
This was attributed to suppression of interference within 



one Dirac valley due to ripples and large defects (Moro- 



zov et al. , 2006 ) . More recently, it was demonstrated that 



both WAL and WL could be observed in graphene un- 



der the proper conditions (Tikhonenko et al., 2009). The 



behavior of the magnetoconductance is shown on Fig. 35 
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FIG. 35 Weak localization and anti-localization in 
graphene. (a) The sample is studied for various carrier den- 
sities (here I, II and III). The QH E show s that the sample is 
indeed two-dimensional (see Sec. VII.D). The sample is an- 
nealed at T 400 K for 2 hours to increase its homogeneity, 
(b), (c) and (d) Variation of the magnetoresistance behavior 
as n and T are changed. A transition occurs between WL and 
WAL and negative magnetoconductance (WAL) is observed 
for low n and high T due to the competition between valley 
scattering and d ephasing processes (see text). Figure from 
( Tikhonenko et a/.t|2QQ9D . 



The data shows WAL (negative magnetoconductance) 
for both decreasing carrier density and increasing temper- 
ature. The increasing temperature reduces the dephasing 
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time of the phase of the electron due to thermal fluc- 
tuations, while decreasing carrier density increases the 
intervalley and intravalley elastic scattering times and 
deflned in the theory of reference ( McCann et al. , 
[2006 ). The latter are roughly temperature independent. 
Thus it is observed that the dephasing time is not the 
only parameter controlling weak localization behavior in 
graphene. Instead, one must consider the ratios = ^ 
and = When they are small, anti-localization oc- 
curs while when they are big, localization occurs. Note 
that there exist combinations of these parameters for 
which the correction to the magnetoconductance is sup- 
pressed. This is shown in Fig. [36] 
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FIG. 36 Phase diagram for weak (ant i-) localization in 
graphene. (a) A schematic representation of interference be- 
tween electron paths. For graphene it is expected to be nat- 
urally destructive and a magnetic field is expected to partly 
restore constructive interference (see text), (b) A phase di- 
agram for weak localization behavior in graphene. It is seen 
that the ratios of the dephasing time to the two valley scatter- 
ing times determine the magnetoresistance behavior. Figure 
from dTikhonenko et a/.|[2009( ). 



In graphene, the quantum corrections to the conduc- 
tance survive at much higher temperatures than for 
2DEG semiconductor structures because the electron- 
phonon scattering is expected to be weak in the system. 
Indeed, weak (anti-)localization can be observed up to 
much high temperatures in graphene, in particular in 



large scale CVD grown graphene (Whiteway et al. , 2010). 



It turns out that electron-electron scattering could be re- 
sponsible for the disappearance of magnetoconductance 



at high temperatures ( Tikhonenko et al. 2009 ) . 



D. Measurements of the quantum Hall effect in graphene 

Here we review the principal features of the QHE in 
graphene and compare them to the conventional results. 
We first discuss the general features of the QHE. 



1. Shubnikov-de Haas oscillations 

The existence of LLs implies the appearance of the so- 
called Shubnikov-de Haas oscillations (SdHO) in the QH 
system. The SdHO consist of a strong modulation of the 
longitudinal resistivity function of the carrier 

density n. The minima of these peaks correspond to a 
completely vanishing longitudinal resistivity. This occurs 
when the chemical potentials of both leads sit between 

Then all the k carrying 
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the LLs, as illustrated in Fig. 
states of the LLs are completely filled so that the elec- 
trons cannot carry current in a LL. The only states that 
carry current are the edge states at the Fermi level (see 
Fig. 37). Furthermore, the magnetic field is such that 



each edge carries current in different directions. Hence, 
the only way for backscattering to occur is for an elec- 
tron to go from one edge to another, which is not feasible. 
Thus scattering is suppressed and no voltage can build 
up, hence pxx vanishes. As the Fermi level is raised, a 
LL may sit at the chemical potential of the leads, which 
allows for some backscattering, hence producing a maxi- 
mum in resistivity. 



Each level provides in 
principle one quantunn of 
conductance eYh 
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FIG. 37 Quantization of Hall conductance in 
graphene. (a) The density of states is peaked at the LLs 
and broadened by the presence of disorder, (b) Edge states 
topologically separated by the magnetic field conduct the cur- 
rent between the leads at fixed potential. In the absence of 
scattering, the quantum of conductance /h carried by each 
channel translates into the quantization of the Hall conduc- 
tivity (Txy (see text). 



2. Conductance quantization 

Associated with the SdHO is the celebrated phe- 
nomenon of Hall conductance quantization. It is ob- 
served that at the values of n for which the zeros in p^x 
occur, the Hall resistivity pxy is constant and forms a 
plateau that is a multiple of h/e^. The integer v is 
said to be the filling factor of the QH state and repre- 
sents the total number of levels below the Fermi energy. 
In fact, there are as many conducting edge channels as 
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there are filled levels, and one can show that each such 
cha nnel should c ontribute one quantum of conductance 



groups (Novoselov et al. 2005 Zhang et al, 2005) and 



( Datta 



1997). In the absence of backscattering, one 
expects each edge to be at the same potential as the lead 
that provides its carriers: this implies that the transverse 
Hall voltage is the same as the potential difference be- 
tween the leads. Therefore the Hall conductance should 
be the same as the sum of the conductance of each of 
these u channels. Thus: 



G 



xy 



(15) 



This phenomenon is known as the quantum Hall effect. 



The reader is referred to (Datta, 1997) for a thorough the- 



oretical treatment of ballistic transport in micron scale 
samples. 



3. Integer quantum Hall effect 

a. Monolayer graphene For the plateaus of the conven- 
tional QHE in semiconductor heterostructures, u is a 
positive integer and a finite number of charge carriers is 
required to occupy the lowest LL. At low magnetic field, 
only the even filling factors appear because the LLs can 
be filled with both spin up and spin down electrons. We 
have the sequence: 



2,4,6,8,10... 



(16) 



Therefore the sequence of i^'s provide information on 
the degeneracies of the electron states or equivalently, on 
their symmetries. If the magnetic field is increased, the 
degeneracies can be completely lifted via Zeeman inter- 
action of the spin with the magnetic field, revealing the 
entire sequence of filling factors. 
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FIG. 38 Typical measurement of the integer quantum 
Hall effect. The Hall resistivity exhibits plateaus quantized 
in exact multiples of Image created by D.R. Leadley, 
Warwick University (1997). 



As a two-dimensional system, graphene is expected to 
exhibit the QHE. It was observed in 2005 by two research 



showed features that were characteristic of the linear dis- 
persion of graphene. A measurement of the conventional 
QHE in graphene is shown in Fig. l39) 




FIG. 39 Integer quantum Hall effect in graphene. The 

carrier density n can be tuned across the neutral point n = 0. 
The absence of a plateau at n = indicates the presence of 
a LL at the neutr al point. T he inset shows bilayer graph ene, 
discussed in Sec. VILD.3.b Figure from ( [Novoselov et al\ 



2005). 



The sequence of observed plateaus is very different 
from that of the conventional QHE. The precise sequence 
is: 



V = ±2, ±6, ±10... 



(17) 



We notice that the carrier density (i.e. the Fermi level) 
of graphene can be tuned through both the valence and 
conduction band of graphene, resulting in negative v for 
holes and positive u for electrons. From the sequence 
we immediately see that each plateau contributes Ae^ /h 
to the conductance, instead of 2e^ /h. This is indicative 
of the expected extra valley degeneracy of the states, in 
addition to the spin degeneracy. Another interesting fea- 
ture of this sequence is that there is no plateau at zero 
carrier density. Thus there is a LL = at the neutral 
point with degenerate holes and electrons as seen from 



Eq. (13). The structure of electron and hole states is also 
completely symmetric. 

Fig. |40] shows a complete lift of the degeneracy of the 
A^ = LL at ultrahigh magnetic fields, confirming explic- 
itly the 4-fold degeneracy of the level ( Zhang et al. , 2006 ) . 
The Zeeman coupling in graphene is too weak to cause 
the lifting of the spin degeneracy and does not explain 



the breaking of the valley pseudospin symmetry ( Zhang 



et al. 2006, 2005). Instead, the lifting of the degeneracies 
is suspected to be caused by enhancement of Coulomb 
and exchange interactions in the sample. This is due to 
the fact that as the field is increased the cyclotron or- 
bits become smaller and the carriers come closer to each 
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other dDean et al] |2QlQa| |Ezawa[ [20081) . integer fill- 



ing fractions can be observed for the cleanest samples 
dDean et a/.l |2010a|). 
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FIG. 40 Complete lifting of the degeneracy of the low- 
est LL with fields from 9 T (circle) up to 45 T (star). The 
expected 4-fold degeneracy is shown explicitly by the presence 
of plateaus at = and = ±1. The plateaus at zy = ±4 
show partial degeneracy lifting of the levels N — ^1. All de- 
generacies can be lift edj n clean samples, see Fig. |44[ Figure 
from dZhang et a/.H20^ . 

Just like the conductance quantization, the SdHO os- 
cillations show unusual behavior. We still observe their 
minima at the positions of the plateaus, as seen in Fig. 39 



A (X 



sinh^^^^P^ 

heB 



(18) 



where rric is the cyclotron mass of the electrons at the 



Fermi level (Novoselov et al. 2005). 



A fit to the data shows that the rapid decay of the 
SdHO implies that rric oc y^, as shown in Fig. 41 Unlike 



in semiconductor devices, the cyclotron mass in graphene 
varies with density. It is possible to show that this depen- 
dence of rUc on n implies the linear dispersion of graphene 
E{k) = hypk and thus the form of the LLs given in 
Eq. ([13| dNovoselov et a/.[|2005[ ). 



b. Bilayer graphene 



The quantum Hall effect can also be 

as 



observed in multilayer graphene. The inset of Fig. 39 
well as Fig. [42] shows measurements carried out on the 
graphene bilayer dNovoselov et al. , 2006). The observed 
sequence is as follows: 



zy = ±4, ±8, ±12, ±16... 



(19) 
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FIG. 41 Determination of the dispersion of graphene 
through SdHO oscillations, (a) The SdHO decay more 
rapidly with temperature as the carrier density is increased. 
The data can be fitted to yield the cyclotron mass rric. (b) The 
cyclotron mass varies with carrier density as ^/n, implying the 
linear dispersion of graphene. Figure from ( ,Novoselov et al.^ 
20051. 



However, their amplitude decays more rapidly with 
temperature as the density n is increased. Theoretical 
studies have shown that for a given n, the amplitude A 



of the SdHO should follow (Castro Neto et al. 2009): 



FIG. 42 Quantum Hall effect in bilayer graphene. As in 

the case of the monolayer, the quantum Hall effect in graphene 
bilayer shows no plateau at zy = 0, indicating the presence of a 
LL at zero carrier density and the relativistic behavior of the 
bilayer. However, contrarily to the monolayer, the iV = level 
is 8-fold degenerate as the height of the step indicates. This 
is consistent with the theoretical prediction of a "parabolic 
Dirac point" at = 0, i.e. massive Dirac fermions. Figure 
from dNovoselov et a/.[|2006 ). 



The absence of a plateau at = is a common char- 
acteristic of the QHE in graphene and bilayer graphene. 
Once again it indicates the presence of a LL at zero en- 
ergy. However, this level is eight times degenerate as can 
be seen from the 8e^//i jump in the Hall conductivity. 
The other levels keep the 4-fold degeneracy of graphene. 
This is what we expect from the additional layer degree 
of freedom of the electron for a given spin and valley 



(Ezawa, 2008 Novoselov et a/., 2006). The degeneracy 



of the lowest LL can be completely lifted by many-body 



interactions in ultra-high magnetic fields (Zhao et al. 



2010), as shown in Fig. 43 
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FIG. 43 The degeneracy of the lowest LL of bilayer 
graphene is lifted by magnetic fields ranging from 9 

T to 35 T. We clearly distinguish for levels above the neutral 
point and 4 others start to be visible below the neutral points, 
accounting f or the expected 8 -fold degeneracy of this level. 
Figure from (iZhao et a/.| 12010). 




FIG. 44 The FQHE in graphene. Muhiple plateaus at 
rational filling factors show that electron interactions give rise 
to new states. Figure from ( [Dean et al. 2010a| ). 
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4. Fractional quantum Hall effect 



FIG. 45 Schematic SdHO structures corresponding 
to different possible symmetries of the lowest LL in 
graphene. (a) Total symmetry breaking, (b) Either the spin 
or the valley symmetry is broken, (c) Full symmetry of the 
spin and valley degrees of freedo m. The latter is the best fit 
to the observations. Figure from (Dean et al. 2010a). 



The extreme two-dimensional confinement in graphene 
is expected to enhance many-body interactions between 
electrons. Since 1982, we know that these interactions 
can manifest through the so called fractional quantum 
Hall effect (FQHE), consisting of the appearance of Hall 
plateaus at fractional (rational) values of the filling factor 
V dTsui et a/.| |l982D . The efforts to observe the FQHE in 



graphene were without success for a long time because of 
the existence of a competing insulating state induced by 
disorder near the neutral point v = {). The improvements 
in the sample fabrication made its observation possible 
in recent years. 

FQH states can be understood as the realization of the 
integer QHE for weakly interacting quasiparticles called 
composite fermions. In a heuristic picture, an even num- 
ber of magnetic flux vortices bind with an electron to 
form an object with reduced effective charge that is a 



Fig. 44 shows a clean and recent 6-terminal measure- 
ment of the FQHE in graphene made on a hexagonal 



boron-nitride substrate at ultrahigh fields (Dean et al. 



2010a). The filling factors appear in the following se- 



quence: 



The measurements for hole carriers were not performed 
for V < —2. The filling factors are all multiples of |, al- 
though I is missing. Furthermore, a peak in the SdHO 
may indicate the emergence of a plateau at = |. The 
effect is much more robust under increase in tempera- 
ture than its counterpart in semiconducting devices, re- 
vealing the suspected enhancement of electron-electron 
interactions in graphene compared to those occurring in 



fraction of the elementary charge e (Dean et a/., 2010a 2DEGs ( Dean et a/.[|2010a[|Du et a/.[|2009[ ). The absence 



Jain 1989). In graphene, the multiple components of the 
electron wave function are expected to give rise to new in- 
teracting ground states with various spatial arrangements 
of pseudospin (textures) that depend on the symmetries 
of the states. These special states are theoretically pre- 
dicted to occur at filling factors that are multiples of — , 



where m is an odd integer (Ezawa, 2008). The FQHE al- 



lows for a characterization of these symmetries and puts 
hard constraints on theoretical models. Different sym- 



metry scenarios are presented in Fig. 45 



= =b| and the presence of = ±^ is consistent with a 
global SU{A) symmetry of the FQH state in the lowest 
LL dDean et a/.||2010aD . 

Previous measurements have been performed on very 
small (<2 /im) 2-terminal ultraclean samples of sus- 
pended graphene ( |Bolotin et al.\ |2009[ |Du et al.\ [2009 ). 
For a discussion of the properties of suspended graphene 
see Sec. |VI.C.2| Their results are presented in Fig. [46[ 
They also obtain the FQHE in the lowest LL at u = ^ 
and z^ = |. While 2-terminal measurements make the 
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results less precise and harder to interpret (for example, 
the contact resistance makes it very hard to extract in- 
formation from the longitudinal resistance data and may 
produce unusual features at v = they possess the ad- 
vantage of being realizable with lower field magnets in 
standard physics laboratories. 




FIG. 46 High quality suspended graphene samples 

used by Bolotin et al. (left) and Du et al. (right). The 
contacts are less than 2 jim. apart in both cases, making a 4- 
terminal measurement impossible. However, the FQHE can 
be observed clearly at relatively low fields ranging from 5 T 
to 12 T. Figures from ( |Bolotin et a/.| [2009] [Pu et a/. | [2009 ). 



The observation of the anomalous quantum Hall effect, 
including its fractional daughter are one of the brightest 
highlights, which propelled graphene as one of the most 
studied materials in recent years. We now move away 
from magneto-transport to discuss the mechanical prop- 
erties in the next section, which are very important for 
future applications in nanodevices. 



This allows us to characterize the quality factor. Young's 
modulus and built-in tension ( Bunch et a/.[|2QQ7i ). One 
can also extract this information from a separate ex- 
periment using an atomic force microscope (AFM) tip 
(Frank et al. 2007). In addition to their use as mechani- 



cal resonators, graphene sheets are robust (impermeable) 
enough to act as a thin membrane between two dissimilar 



environments (Bunch et al. 2008). 



A. Overview of the harmonic oscillator 

In order to study the mechanical properties of 
graphene and graphite sheets, one needs to comprehend 
the basics of a harmonic oscillator. A classic example is 
a mass attached to a spring. In an ideal situation, the 
spring obeys Hooke's law, but for large displacements, 
the spring does not follow Hooke's law as the system is 
subject to damping, denoted by 7, which dissipates the 
vibrational energy of the system. In addition, a driving 
force F with frequency uj is necessary to drive oscillations 
to the system. Therefore, the equation of motion for a 
damped harmonic oscillator is given by: 



d X dx 9 . . 

m— y + 7—- + tuuOqX — h cos[ujt). 

(JjL (JjL 

The general solution is given by: 



(21) 



VIII. MECHANICAL PROPERTIES 

Microelectromechanical systems (MEMS) have been 
deployed to perform common tasks such as opening and 
closing valves, regulating electric current, or turning mir- 
rors. This is realized by employing microscopic machines 
such as beams, cantilevers, gears and membranes. These 
MEMS are found in many commercially available prod- 
ucts from accelerometers found in air bag deployment 
systems, gyroscopes in car electronics for stability con- 



trol and ink jet printer nozzles ( Craighead] 2000| ). The 
nanoscopic version of MEMS, nanoelectromechanical sys- 
tems (NEMS), demonstrate their own advantages in en- 
gineering and fundamental science in areas such as mass. 



force and charge detection (Ekinci and Roukes, 2005). 



All these electromechanical devices are functional only 
in response to an external applied force. The utmost 
limit would be a one atom thick resonator. Robustness, 
stiffness and stability is thus important when reaching 



this limit (Bunch et al. 2007) 



The fabrication of graphene-based mechanical res- 
onators is still under active development. Presented here 
is an overview of experimental results for single and mul- 
tilayer graphene sheets placed over predefined trenches. 
A simple drive and detection system is used to probe 
the mechanical properties of these graphene resonators 
in order to extract the fundamental resonance frequency. 



X = Acos{ujt — 0) 
where is the phase shift, the amplitude 



A = 



Q 



muj. 



41 ^^i^J 



(22) 



(23) 



and the quality factor 



Q 



mujo 



(24) 



For a given peak with resonance frequency uq = 27r/o, 
the peak amplitude is QF/mujQ and its full- width half 
maximum (FWHM) determines the Q factor, given by 

fo/Q- 



B. Tuning resonance frequency by electrical actuation 

The mechanical resonators are exfoliated graphene 
sheets that are suspended over trenches on Si02/Si sub- 
strates. These micron-sized sheets are doubly-clamped 
beams which are secured to the Si02 surface via van der 
Waals forces. A gold electrode defined by photolithogra- 
phy is used to electrically actuate the graphene resonator. 
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FIG. 47 Schematic of a suspended graphene resonator 



for electrical actuation (Bunch et al. 2007) 



filling the role of the driving force F from Eq. (21). All 
resonator measurements are performed in ultra-high vac- 
uum at room temperature. The actuation is an applied 
alternating electric field used to drive the resonant mo- 
tion of the beam. A capacitor is then formed between 
the bottom gate electrode and the contacted graphene, 
shown in Fig. 47 A voltage Vg applied to the capacitor 
induces electric charges onto the beam. A small time- 
varying radio frequency (RF) voltage 6Vg at frequency 
/ is used to electrically modulate the beam on top of a 



constant DC voltage (Bunch et al 2007): 



V,=Vf^ + SV,. (25) 
The resulting electrostatic force is then given by: 



mode remains unaffected. A possible explanation for the 
tuning may originate from an electrostatic attraction to 
the gate which increases the tension from stretching the 
graphene. 
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FIG. 48 Resonance spectrum taken with an electro- 
static drive. Plot of the amplitude versus frequency of the 
(a) fundamental mode at 10 MHz and (b) the higher mode at 
35 MHz with increasine; 6Vg while Vg^^ remains fixed. Inset: 
Plot of the resonant peak amplitude with increasing SVg. De- 
pendence of the (c) amplitude and (d) frequency as function 
of Vf^ at fixed 6Vg where the solid squares and triangles rep- 
resent the fundamental and higher mode respectively ( |Bunch| 
et al.\\2007\ . 



2 dz ^ 



where z is the distance between the graphene and the 
gate electrode. 

The response due to a varying RF voltage is monitored 
by a 632 nm He-Ne laser focused on the resonator. As 
a result, the light creates an interference between the 
suspended graphene sheets and the silicon back plane. 
Variations of the reflected light intensity are monitored 
by a photodiode. 



Plugging Eq. (26) into Eq. (23), one can speculate on 



the behavior of the resonant peak as a function of applied 
voltage. The data for a few-layer graphene stack is shown 



in Fig. 48 for electrical drive on resonance. One can see 
that the amplitude and frequency increase linearly with 
5Vg at a flxed for the higher mode (Fig. [48p) while 
the frequency of the fundamental mode (Fig. |48|A.) does 
not change signiflcantly. In addition, if Vg is increased 
while 8Vg is flxed (Fig. 48 3), both the amplitude of the 



fundamental mode and that of the higher mode increase 



linearly as expected from Eq. (23). As for the frequency 
(Fig. pi 



it seems that there is evidence of positive 
tuning for the higher mode but that the fundamental 



C. Resonance spectrum by optical actuation 

Another technique to actuate vibrations in resonators 
is an optical drive. A diode laser is shined on the 
suspended graphene in which the graphene actuates 
by itself via thermal expansion and contraction. This 
motion is controlled by the intensity of the diode laser 
set to a frequency /. This is then monitored by a 
photodiode using a He-Ne laser through the procedure 



given in the previous section (Bunch et al. 2007). 



Fig. [49JA. shows the resonance spectrum of a few-layer 
graphene sheet suspended over a Si02 trench. Several 
resonant peaks were measured; the flrst peak is the 
fundamental vibrational mode which will be the main 
focus of the discussion. The fundamental peak is the 
most pronounced peak, so it is easy to obtain the 
mechanical characteristics of the resonator. Analyzing 
higher modes is unnecessary as they depend on the 
fundamental frequency /q. The most appropriate flt to 
extract the resonance frequency /o is a Lorentzian flt, 
giving /o to be 42 MHz with a quality factor Q = 210. 
Fig. [49^ shows similar results for a single layer graphene. 
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FIG. 49 Plots showing the resonance spectrum of 
graphene and graphite resonators, (a) Resonance spec- 
trum showing the fundamental mode and higher modes for 
a 15 nm thick multilayer graphene sheet taken with optical 
drive. Inset: Optical image of the resonator. Scale bar = 5 
/xm. (b) Amplitude vs. frequency for a single layer gr aphene. 
The red curve represents the Lorentzian fit of the data (Bunch' 
et a/.||2007|). 



/o = 70.5 MHz and Q = 78. The same measurements 
were repeated for 33 resonators of different geometries 
with thicknesses ranging from 1 atomic layer to 75 nm. 
They are plotted in Fig. [50] The data shows that the 
fundamental frequency /o largely varies from 1 MHz 
to 166 MHz with quality factors ranging between 20 
and 850. Knowing the fundamental resonance /o, the 
following equation for a doubly-clamped geometry is 
used as a reference to extract the Young's modulus of 



the sheets (Weaver et al. 1990): 



\ 




Q.blA^S 
pLH ' 



(27) 



where E is Young's modulus, S is the tension per 
width, p is the mass density, t and L are the thickness and 
length of the suspended graphene sheet, and ^4 is a ge- 
ometrical constant; A = 1.03 for doubly-clamped beams 
and 0.162 for cantilevers. The parameter S represents the 
built-in tension on the graphene sheet which may come 
from the fabrication process or from the van der Waals 
interaction between the substrate and graphene. Assum- 



ing that the tension is sufficiently small, Eq. (27) pre- 
dicts that the fundamental frequency /o scales as t/L'^. 
Fig. [SOj A is a plot of the resonance frequency for res- 
onators with t > 7 nm (solid squares) and t < 7 nm (hol- 
low squares) as a function oit/L'^ and the slope gives the 
Young's modulus E. Taken the density for bulk graphite 
as p = 2200 kg/m^, the Young's modulus of the graphene 
sheets plotted as dashed lines correspond to values be- 
tween 0.5 and 2 TPa. These values correspond well to 
the Young's modulus of 1 TPa of bulk graphite ( |Kelly[ 
[l981j). This is the highest modulus resonator to date. 
This is in stark contrast to the 12-300 nm thick Si can- 
tilevers (Sazonova et al. 2004 ) which achieve values rang- 



triangles in Fig.|50j^. For graphene resonators with t <7 
nm, the data points are not as linear as for the thick res- 
onators. To elucidate these results, Scharfenberg et al. 
studied the elasticity effect of thin and thick graphene 
samples on a corrugated elastic substrate such as poly- 
dimethylsiloxane (PDMS). For thin graphene samples, it 
adheres fully to the substrate so that the graphene follows 
the topography of the substrate. As for thicker samples, 
it flattens the corrugated substrate. These observations 
suggest that thin graphene flakes, despite their size, are 
highly sensitive to tension, making them easier to deform 



than the thicker samples (Scharfenberg et al. 2011). 
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FIG. 50 Measurements of the Young's modulus and Q 
factor for doubly-clamped beams, (a) Plot showing the 
frequency of the fundamental mode of all the doubly-clamped 
beams and Si cantilevers versus t/L^. The cantilevers are 
shown as solid triangles. The doubly clamped beams with 
t > 7 nm are indicated as solid squares while those with 
t < 7 nm are shown as hollow squares. The solid line is the 
theoretical prediction with no tension, (b) Plot of the quality 



factor as a function of the thickness for all resonators ( Bunch 



et al. 2007). 



The quality factor is an important parameter to con- 
sider when dealing with resonators. It gives a good indi- 
cation of the resonator's sensitivity to external perturba- 
tions. A high Q is thus essential for most applications. A 
plot of Q as a function of the thickness for all graphene 
resonators is shown in Fig. [50^. It seems that there is 
no clear dependence of quality factor on resonator thick- 
ness observed and the Q factor is on the order of hun- 
dreds. These Q factors are lower than diamond NEMS 

(Q^ 



^2500-3000) dSekaric et al.\ |2002[ ) and significantly 
lower than high tensile Si3N4 (Q:^200,000) ( |Verbridge 



et a/.[[2QQ6 ). Thicker variants of graphene, like graphene 



oxide (Q:^4000) ( [Robinson et al.\ |2008b|) and multilayer 
epitaxially grown graphene (Q^^IOOO) ( [Shivaraman et oT 



2009) have shown higher quality factors. Despite the us- 
ability of graphene as a mechanical resonator, the source 
of dissipation is still poorly understood. Recent progress 
on graphene resonators suspended over circular holes on 
SiNx membranes have shown some improvements on the 



Q factor. Thus, Fig. [51] demonstrates a clear dependence 
of quality factor on the resonator diameter. The highest 
Q factor is about 2400 for a resonator with a diameter of 



ing between 53 - 170 GPa, and which are plotted as solid 22.5 am (Barton et al. , 2011). 
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FIG. 51 Resonance frequency and Q factor of 
graphene membranes, (a) Fundamental frequency / as 
a function of diameter D for 29 graphene membranes. The 
red hne is a fit of the data in which / - jj-o.9±o.i^ ^^^^ 
Quality factor of membranes as a function of diameter. Inset: 
the highest Q factor peak observed fitted with a Lorentzian 
indicating a value of 2400 ( ^Barton et al.\ 2011). 



D. Measurements by atomic force microscopy 

In the previous subsections, the Young's modulus is 
measured by considering the tension being small, neglect- 



ing the second term in Eq. (27) for a doubly-clamped 



beam. The built-in tension in the suspended graphene 
may indeed play a role in the mechanical measurements, 
a plausible explanation to the scattered data found in 
Fig. 50 for thin graphene resonators. To measure the 
built-in tension, an AFM tip with a calibrated spring 
constant is pressed onto the suspended graphene (Frank 
et al 2007). The spring constants on the sheet are 
then extracted to calculate the Young's modulus and 
built-in tension. With an applied static force for a 
doubly-clamped beam and employing the relation /o = 



(l/27r) Y^/cg/m and Eq. (27), the resulting effective spring 
constant is: 



ks = 16.23EW 



L 



4.93 



L' 



(28) 



where w is the width of the sheet. 
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FIG. 52 Schematic of a deflected AFM tip v^hile push- 
ing on the suspended graphene. r]graphene is the height 
that the tip is deflected by the graphene and Azgraphene is 
the height pushed by the tip (Frank et al.\ 2007). 



Fig. [52] shows the schematic of an AFM tip pushing 
down the suspended graphene sheet. The deflection of 



the tip then helps to find the effective spring constant of 
the suspended sheet. The chosen AFM tip has a spring 
constant 2 N/m which enables the graphene to be de- 
flected by a detectable amount. The tip is pushed slowly 
against the sheet and the relation between tip displace- 
ment and the position of the piezo is then plotted in 
Fig. 53 As the AFM tip comes into contact with the 



suspended graphene, the cantilever is pulled down onto 
the surface resulting in a dip in the deflection as observed 



in Fig. 53 From the tip's spring constant, a graph of the 
force exerted on the tip versus the displacement of the 
graphene sheets can be extracted as plotted in Fig. [53) 3. 
The displacement can be written as: 



— Vgraphene ^^graphe 



(29) 



where rjgraphene IS the deflection of the tip, Zpie^o is 
the location of the piezo moving the tip, and Azgraphene 
is the deflection of the graphene sheet. From Hooke's 
law, the slope yields the effective spring constant of the 
suspended graphene sheet. 
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FIG. 53 Measurement of the graphene's spring con- 
stant, (a) Curve obtained from the deflection of the AFM 
tip by pushing down the suspended graphene sheet. The right 
axis represents the force corresponding to the tip displace- 
ment, (b) Plot of the force as a function of the displacement 
of the graphene sheet ( Frank et al. 2007 ) . 



Fig. [54] shows the plot of the spring constant as a 
function of the dimensions of the suspended graphene 
sheet for eight samples. The spring constant varies 
between 1-5 N/m. A linear fit is then used to extract the 
built-in tension and Young's modulus of the suspended 



graphene sheets. From Eq. (28), the slope suggests that 
a Young's modulus E of 0.5 TPa, in good agreement 
with suspended graphene actuated optically and bulk 
graphite. The offset of the linear fit gives a tension 
of 300 nN, suggesting that the built-in tension of all 
suspended graphene sheets are on the order of nN. 

Other work on AFM nanoindentation was realized on 
graphene suspended over micron-sized circular wells to 
investigate the elasticity properties of graphene. Indent- 
ing defect-free graphene with an AFM tip similar to 



Fig. 52, one can probe the elastic stress-strain response. 
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FIG. 54 Plot of the spring constant in the center of 
the suspended region versus w{t/L)^ for eight different 
samples with various thicknesses. The linear fit provides 
information about the built-in tension and Young's modulus 
of the graphene sheets (Frank et al. 2007). 



A non-linear elastic response model is fitted to find that 
the stiffness is about 340 N/m with intrinsic breaking 
strength of 42 N/m ( |Lee et al\ |2QQ8 ). 



As discussed in this section, graphene as an ultra-thin 
membrane with extremely high bulk modulus is very 
promising for potential applications in mechanical sys- 
tems, if the source of dissipation can be reduced sub- 
stantially, which should lead to a significant increase in 
Q factor. We now turn to applications of graphene as 
transistors. 



duced, more and more issues (commonly known as short 
channel effects) begin to appear ( [Frank et a/. , 2001). 
These issues include but are not limited to such prob- 
lems as hot electron effects, velocity saturation effects, 
and punchthrough effects. It has been suggested that 
graphene, due to its monoatomic thickness, would reduce 



these parasitic effects (Schwierz 2010). Graphene could 



be included in the channel region of transistors and thus 
provide a high mobility channel which would help to re- 
duce the effect of short channel effects. One problem is 
that graphene does not possess a bandgap since the con- 
duction and valence bands touch at the Dirac point. This 
reduces the on-off ratio of the transistor by several orders 
of magnitude, such that a graphene FET will typically 
have an on-off ratio < 10. This ratio is unacceptable 
if one wants to replace CMOS technologies where ratios 
> 10^ are commonly achieved. It has been proposed 
by the International Technology Roadmap of Semicon- 
ductors (ITRS) that a ratio of 10^ would be required 
for logic applications. Such a ratio could be achieved in 
graphene only by opening a bandgap at the Dirac point in 
order to suppress the band-to-band tunneling ( [Fiori and] 
lannaccone 2009). A bandgap of 0.4 eV could achieve 



the desired on-off ratio (Schwierz 2010i). Several tech- 



niques are being envisioned for engineering a graphene 
bandgap. The following sections provide an overview of 
these methods. 



1. Bilayer graphene 



IX. GRAPHENE TRANSISTORS 

This section provides an overview of current experi- 
mental work being pursued toward the development of 
graphene field effect transistors (GFET). The two most 
common varieties of transistor are the logic transistor 
and the analog transistor. The former is characterized, 
among other things, by a high lon/^off ratio to ensure 
low energy consumption (of the off state) and to main- 
tain a high logic interpretation yield. The latter is char- 
acterized primarily by its cutoff frequency, and is often 
used in high frequency applications as an amplifier. Both 
types of devices, as we shall see, could benefit from the 
remarkable electrical properties of graphene. 



A. Logic transistors 

For more than forty years CMOS technology has dom- 
inated the logic transistor industry with the fabrication 
of MOSFETs. A general trend over the years has been 
to reduce the length of the transistor's gate in order to 
achieve higher transistor densities but also increased per- 
formance produced by the higher electrical fields in the 
channel region. However, as the size of the device is re- 



The first technique that we shall investigate is the use 
of bilayer graphene. Theoretical calculations have pre- 
dicted the possibility of a bandgap opening in this con- 
figuration ( [McCann] 2006). Such a structure allows the 
opening of a bandgap by breaking the symmetry of the 
bilayer stack (Bernal stacking) with the application of a 
transverse electrical field. One of the first investigations 



of this approach is given in reference ( Ohta et al. , 2006 ) . 



The authors used bilayer graphene on a SiC substrate. 
In such a structure, light doping is provided to the bot- 
tom graphene layer by the substrate and further doping 
can be introduced to the top graphene layer by the ad- 
sorption of potassium atoms. This setup was designed 
and fabricated to facilitate band structure analysis using 



the ARPES technique. Fig. [55] presents the experimental 
band structure for increasing amounts of potassium dop- 
ing and a comparison to theoretical calculations. One 
may notice the good agreement between probed and the- 
oretical bandstructure. A maximal bandgap of approxi- 
mately 0.2 eV is observed using this technique. 

The results presented above show the potential of bi- 
layer graphene for bandgap engineering. Several groups 



used this material for designing FET devices (Oostinga 



et 



2008). In these devices doping is induced using 



the field effect instead of using potassium atoms. Fig. 56 
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FIG. 55 Bandgap evolution as a function of the potas- 
sium doping. The doping level is controlled using the potas- 
sium adsorption of bilayer graphene. Shown from left to right 
is a diagram of the evolution of the bandgap for an increasing 
doping level per unit cell ( |Ohta et a/.[|20Q6[ ). 



presents a schematic representation of such a device. The 
double gate configuration allows one to control both the 
Fermi position and the bandgap size. 



On one hand, this approach makes it easy to control 
the position of the transistor's elements (e.g. channel, 
gate, electrons, etc.). On the other hand, making the 
width of a nanoribbon arbitrarily small is limited by 
the lithographic process. Furthermore, the edge quality 
of the patterned nanoribbons has been found to be 
rough, leading to increased scattering and decreased 
performance ( |Han et al. |2007l). Fig. 57 presents a 



graph of the bandgap energy (Eg) as a function of the 
width of the lithographically patterned nanoribbons at 
T = 4.2 K. Again decreased performance is correlated 
with increasing temperature. One can notice a maximal 
bandgap of '^200 meV for a nanoribbon width of ~15 
nm. These results are similar to the ones obtained for 
bilayer graphene. 




FIG. 56 Schematic of a GFET device. Graphene is de 
posited on an isolating material (e.g. Si02) after w hich regu- 
lar photolith ographic patterning can be performed ( Oostinga 
et al. 



2008) 



After characterization of the device, a maximal 
lon/Ioff ratio of approximately 100 could be achieved 
at 4.2 K. Decreasing performance was observed with in- 
creasing temperature. These low on-off ratios are ex- 
plained by the insufficiently large bandgap. Theoreti- 
cal calculations for bilayer graphene predict a maximal 
achievable bandgap that would be below the 0.4 eV gap 
needed for an on-off ration over 10^. These predictions 
are corroborated with experimental results. Therefore, 
bilayer graphene seems to be impractical for logic appli- 
cation. 



2. Graphene nanoribbons 

A second technique used for engineering a bandgap in 
graphene is the use of graphene nanoribbons. The idea 
is to provide additional charge carrier confinement in 
the plane of the graphene sheet, effectively making it a 
one-dimensional structure. This additional confinement 
can lead to the opening of a transport gap. There are 
two common techniques for producing nanoribbons. The 
first technique is the lithographic patterning of graphene. 
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FIG. 57 Energy bandg ap as a function of nanoribbon 
width at T = 4.2 K ( ,Han et a/. [ [2007 J . 



A second approach for producing nanoribbons is chem- 
ical synthesis (Li et al. 2008). This method, unlike 



lithographic patterning, produces smooth edges on the 
nanoribbons, while achieving widths in the sub- 10 nm 
range. However, this method leads to a distribution of 
nanoribbon sizes and makes their positioning difficult. 
Scalability of such a technique is thus hard to achieve. 
Fig. [58^ presents an I-V curve for a 5 nm wide nanorib- 
bon in a GFET device. One may notice that even at room 
temperature lon/Ioff ratios above 10^ are observed. This 
behavior is explained by the size of the bandgap, which 
reaches 0.4 eV, as illustrated in Fig. [58)3. For such small 
widths, carrier confinement is appreciable leading to an 
asymptotic behavior of the bandgap as a function of 
width. 
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FIG. 58 Chemically derived nanoribbons ( Li et al] 



3. Other techniques 

This section presents three other techniques for 
engineering a bandgap in graphene: strained graphene, 
nanomesh graphene and patterned hydrogen adsorption 
graphene. 



a. Strained graphene Several studies suggest that a 
bandgap would be induced in graphene by the applica- 
tion of strain. By keeping track of the G and 2D bands 
while applying strain on a graphene sheet, one can 
monitor the amount of strain induced. Strain of 0.8% 
has been experimentally measured using this technique 



(Ni et al. 2008b). Theoretical calculations predict that 



such an amount of strain would lead to a 300 meV 
bandgap. However, experimental evidence of this gap 
has not been demonstrated. 



lon/Ioff ratio of the order of 100 at room temperature 
for patterned holes of ^7 nm in width. The on-off 
ratio tends to increase with decreasing neck width 
of the holes. A major advantage of this technique 
is its scalability. Work remains to be done in order 
to study whether an on-off ratio of > 10^ can be achieved. 



Graphene grown on Ir give rises to a superperiodic po- 
tential, also called a Moire pattern. Such a pattern re- 
mains after the hydrogen adsorption and it alters the 
electronic properties of graphene leading to a bandgap 
opening. Fig. 59 presents ARPES measurements of the 



hydrogen adsorbed graphene. One can easily see the the 
opening of a large bandgap around the Fermi level. This 
structure is also stable at room temperature. 




FIG. 59 Valence band evolution. Exposition time from 
left to right goes from sec (clean graphene) to 30 sec and 
finally 50 sec. A bandgap of almost 0.5 eV is observed af- 
ter 50 sec exposition a result that is in go od agreement with 
theoretical predictions (Balog et al. 2010). 



4. Techniques comparison 

Presented below is a summary table that presents 
achieved results for the different bandgap engineering 
methods presented in this section. 



Method 



Bandgap Ipn/Ipff Ref. 



Bilayer 



0-0.25 eV -100 



(Ohta et al. 2 0£6||Oost-| 
inga e t aLj |200 



Nanoribbons 0-0.4 eV > lO"" 



( Han et 



b. Nanomesh graphene Another novel way of inducing a 
bandgap in graphene is to pattern holes in a graphene 
sheet (Bai et al.\ 120101). This technique led to an Nanomesh 



(Li et al. 200< 



0^ 



20071 



Strain 



0.3 eV* 



NA 



(Ni et al. 2008b^ 



NA 



10 — 100 (Bai et al. 2010) 



Hydrogen 
adsorption 



> 0.4 eV NA 



( [Balog et"ar||2010| ) 



TABLE IV Bandgap engineering techniques *No exper 

mental data available 



B. Analog transistors 



c. Patterned hydrogen adsorption This novel technique 
makes use of half-hydrogenated graphene on Ir(lll) (Ba- 
log et al, 2010). Theoretical DFT simulation predicts 



that this technique can open a bandgap of 0.43 eV in 
graphene, thus reaching our desired value of 0.4 eV. 



A second type of transistor is the analog transistor. 
These devices do not require a non-conducting off-state 
like logic transistors, and therefore do not require a 
bandgap. They are used in radio frequency applications 
and are characterized by their cutoff frequency ft. The 



37 



high mobihty of graphene possibly leads to very high 
cutoff frequencies, making graphene an attractive mate- 
rial for designing these transistors. One major challenge 
for the integration of graphene in these transistors is the 
preservation of high mobility while integrating graphene 
in a device. 

Since 2008 there has been increased research interest 



in the GFET analog transistor (Meric et al , 2008). Some 



of the initial work showed that a 14.7 GHz transistor was 
achievable using graphene with a 500 nm device length. 
Following these results, a 26 GHz device was fabricated 
using a shorter gate length of 150 nm and mechanically 
exfoliated graphene (Lin et a/., 2009). The cutoff fre- 



quency of GFET devices has proven to follow its well- 
known behavior for this type of transistors given by: 



(30) 



27tCg 



Shorter gate length increases the drift velocity of 
charge carriers by increasing the electrical field. It also 
leads to a reduced distance between the source and drain 
contacts, which also reduces the transit time across the 
FET. It was observed that ft cx 1/Lq , so reducing 
the gate length from 500 nm to 150 nm resulted in 
the cutoff frequency rising from 3 GHz to 26 GHz (Lin 



et al 



2009|). For these devices the effective mobility 

^. This value is 



was estimated to be 400 cm^V ^s 



well below the 10000 cm^V ^s ^ mobility of exfoliated 
graphene at room temperature that could be achieved 



( jNovoselov et al. , 2004). Indeed, the effective mobility 
in graphene is highly sensitive to its environment and 
decreases as the amount of scattering elements and 
defects are increased (e.g. insulating substrate, gate 
oxide, source/drain contacts, etc.). 

A similar device configuration using a higher quality 
graphene (graphene grown on SiC) achieved a higher 
effective mobility of 1000-1500 cm^V~^s~^. For this 
device, a 100 GHz cutoff frequency been measured with 



a gate length of 240 nm (Lin et al 2010). 



Lately, other device configurations that minimize scat- 
tering with the graphene sheet are being engineered. 



Fig. 60 shows a GFET device using a nanowire gate. 

This configuration minimizes scattering and defects in 
the graphene sheet and allows one to achieve cutoff fre- 
quencies in the 100-300 GHz with gate lengths in the 
200 nm range. Further work is still required in order to 
achieve higher effective mobility. In the following section 



we shall see how the optical properties of graphene are 
of a great interest for optoelectronic devices. 

.pt ' , ^ ^ 




FIG. 60 GFET schematics using a nanowire gate (Liao 



et al. 



2010). 



X. OPTOELECTRONICS 



A. Optical properties 



The properties of graphene make it an attractive 
choice for use in optoelectronic devices. In particular, 
graphene's high transparency, low reflectance, high car- 
rier mobility and near-ballistic transport at room tem- 
perature make it a promising choice for transparent elec- 
trodes. Optically, single-layer graphene (SLG) has an 
unusually high absorption (given its thickness) that can 
be described in terms of the fine structure constant, a 
(a fundamental physical constant that describes the in- 
teraction of matter and electromagnetic fields). Light 
transmittance T through free-standing graphene can be 
derived using the Fresnel equations for a thin film with 
a universal optical conductance of Go = e^/4fi: 



T = (1 + 0.57ra)"^ ^ 1 - Tra ^ 0.977 



(31) 



where a = e^/(47reo^c) = Go/(7reoc). For the most 
part, each SLG will contribute an absorbance of A = 1 — 
T ^ 7ra ~ 2.3% to visible light. Because graphene sheets 
behave as a 2-dimensional electron gas, they are optically 
almost non-interacting in superposition (though the same 



cannot be said of their electronic interaction (Luican 
\et al.\ |2011 )), so the absorbance of few-layer graphene 
(FLG) sheets is roughly proportional to the number of 
layers. The absorption spectrum of graphene from the 
ultraviolet to infrared is notably constant around 2-3% 



absorption, as shown in Fig. [61fi , compared to some other 
materials. The reflectance is very low at 0.1%, though 
increases to 2% for 10 layers dCasiraghi et al. 2007). 



Graphene has found its way into many optoelectronic 
and photonic devices, a portion of which will be cov- 
ered here. Other uses, including photodetectors, touch 
screens, smart windows, saturable absorbers and optical 
limiters have recently been covered in detail elsewhere 



(Bonaccorso et al. 2010). 
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FIG. 61 Exploring the utility of graphene for thin conducting electrodes. Reprinted from (Bonaccorso et al. 2010). 

Theoretical ranges are shown bound by black lines, as calculated using typical values with Eq. (32). (a) Transmittance as a 



function of wavelength for graphene, ITO and two other metal oxides, as well as SWNTs. (b) Rs values as a function of thickness, 
illustrating the potential of graphene to achieve substantially lower Rg for a given thickness, (c) Transmittance as a function 
of Rs for the same materials, (d) Transmittance as a function of Rs for a number of graphene fabrication methods: triangles, 
CVD; blue rhombuses, micromechanical cleavage (MC); red rhombuses, organic synthesis from polyaromatic hydrocarbons 
(PAHs); dots, liquid-phase exfoliation (LPE) of pristine graphene; and stars, reduced graphene oxide (rGO). 



B. Transparent conducting electrodes 

Transparent conducting electrodes (TCEs) are criti- 
cal components in a wide variety of devices, from those 
found in academic and industrial settings to commer- 
cial devices that define modern technology. They func- 
tion primarily by injecting or collecting charge depend- 
ing on the purpose of the device, and are highly trans- 
parent across some part of the electromagnetic spectrum 
(most often to visible light). The most widely used ma- 
terial for TCEs is indium tin oxide (ITO), an n-type, 
wide-bandgap semiconductor {Eg = 3.75 eV), which typ- 
ically consists of 90% indium(III) oxide doped with 10% 
tin(IV) oxide by weight. The tin atoms function as n-type 
donors. ITO has favorable electronic and optical proper- 
ties for most TCE applications, but suffers from practical 
limitations. It is relatively inflexible and fragile, and a 



limited supply of indium makes it increasingly costly for 
large-scale production. Though it is highly transparent 
across most of the visible range and near-infrared, it be- 
comes increasingly opaque at UV and reflects IR wave- 
lengths, because of band-to-band absorption (excitation 
of an electron from the valence to the conduction band) 
and free carrier absorption (excitation within the conduc- 
tion band), respectively (Keshmiri et al, 2002). Despite 



the utility of ITO in optoelectronic devices, alternatives 
are sought due to the aforementioned reasons. One pop- 
ular inorganic substitute is doped zinc oxide, which has 
comparable properties, though suffers from limited etch- 
ing capability due to acid sensitivity, among other things. 
It is typically doped with aluminum, gallium or indium. 
The list of inorganic semiconductor materials is exten- 
sive, but ITO has become the standard to beat for many 
applications. 
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In contrast to inorganic materials, organic TCEs are 
characterized by low cost, flexibility and stability, though 
do not achieve the same level of charge mobility as 
inorganic materials. Many varieties of intrinsic con- 
ducting polymers (ICPs) have found use in optoelec- 
tronics. In particular, their electroluminescence, flex- 
ibility and mechanical strength makes them particu- 
larly suitable for OLEDs and touch screens, though 
processability is a known issue. PEDOT (poly(3,4- 
ethylenedioxythiophene)) is a conducting polymer that 
is widely used as a transparent conductive film in the 
form of PEDOTipolystyrene sulfonic acid (PSS) disper- 
sions (due to the insolubility of PEDOT in water). Like 
other ICPs, the conductive properties of PEDOT arise 
from a chain of conjugated double bonds, and by con- 
trolling the extent of 7r-electron cloud overlap, the PE- 
DOT band gap can be varied from 1.4 to 2.5 eV in a 
manner analogous to doping (earning such polymers the 
nickname "organic metals") ( [Groenendaal et aT , 2000). 
For organic semiconducting materials, the highest occu- 
pied molecular orbital (HOMO) and lowest unoccupied 
molecular orbital (LUMO) are analogous to the valence 
and conduction band of inorganic semiconductors. 



C. Graphene as TCE 

Two of the most vital parameters for TCE materials 
are low sheet resistance Rs and high optical transmit- 
tance, T. Sheet resistance (in units of Q/n) is directly 
obtained using a four-terminal sensing measurement (see 
Sec. VI), and can be thought of as resistance per as- 



pect ratio. Under most circumstances, graphene matches 
or exceeds the transmittance of competing materials, 
though the sheet resistance has proven somewhat less 
predictable, and varies depending on production method, 
with CVD being the closest to optimal for optoelectron- 
ics (highest T and lowest Rg). This is illustrated graph- 
ically in Fig. [61] T and Rg in doped FLG are related 
fairly simply by: 

T=l + ^— (32) 

ZUs (J2D 

with optical conductance Go as defined above, and 
where Zq = I/cqc = 377 Q is the free-space impedance, 
eo is the permittivity of free space and c is the speed of 
light. Note that this relation is not directly related to the 
number of layers, though the 2-dimensional conductivity 
is given by a2D = ^Me^, where n is the number of charge 
carriers. By taking a range of realistic values for FLG, it 
is possible to determine a range of T and Rs values, as 
shown in Fig. 61 For CVD graphene, n can vary from 
10i2_ioi3 and /ie from 1000-20000 cm^V-^s-^ No- 

tably, by increasing the thickness, and thus the number 
of charge carriers, it is possible to achieve a lower Rg 
value at the expense of transparency. 



1. Work function, Schottky junctions and photodetectors 

An important factor in the design of optoelectronic de- 
vices is the work function of graphene (^g)- A difference 
in the work function between materials at an interface 
results in charge transfer. A particular case of this is 
a metal/semiconductor interface, referred to as a Schot- 
tky junction, where unmatched work functions result in 
a potential barrier with rectifying characteristics. If the 
rectifying effect is absent (depending on the metal work 
function and band gap of the semiconductor), it is re- 
ferred to as an ohmic contact. Schottky diodes differ 
from p-n junction diodes by lower forward voltage and 
faster switching speeds, allowing for high frequency sig- 
nals. 

Despite the absence of a band gap in graphene, 
it can function as the semiconductor active element 
for a metal/graphene photodetector (|Mueller et al. , 
2010). Normally, photoexcitation of graphene results 



in electron- hole pairs that quickly recombine. Near the 
Schottky junction, photogenerated charge is separated 
by the electric field of the Schottky barrier and a pho- 
tocurrent results. In graphene, both electrons and holes 
have high mobility, providing an advantage over conven- 
tional semiconductors. In a preliminary test of an inter- 
digitated metal-graphene-metal photodetector, no signal 
degradation was observed up to the frequency limit of 
the measurement system (40 GHz), and it is believed 
to be capable of up to 0.6 THz, ultimately limited by its 
RC constant ( Avourisj 2010[ ). Using multiple metal types 
(with high and low work function) produces different dop- 
ing, and by incorporating both, the photodetection area 
is maximized. This device was able to reliably detect 
optical data streams of 1.55 jiia light at 10 GBits/s. 



2. Light-emitting diodes 

Light-emitting diodes (LEDs) function by radiative re- 
combination (also known as the Lossev effect), wherein 
recombination of electron-hole pairs can result in the 
emission of photons. They offer a number of advantages 
over incandescent light sources including lower energy 
consumption and longer lifetimes. The general construc- 
tion of an LED includes an electroluminescent semicon- 
ductor material (organic or inorganic) that emits light 
upon recombination and decay of charge injected from 
surrounding electrodes. The most commonly used an- 
ode material is ITO. In addition to the limitations stated 
above, it and other inorganic components may degrade 
over time, leaking ions into the active layer and compro- 
mising the proper function of the device. 

A number of recent publications report gallium ni- 
tride (GaN)-based UV LED devices that use FLG as a 
transparent electrode (Jo et al., 2010a[[Kim et al. 2010). 
Graphene has better transparency to UV/blue light than 
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ITO, and its thermal conductivity (particularly for a 
given thickness) is advantageous for improved heat distri- 
bution and extension of the lifetime of the device. Nitride 
(including gallium) LEDs are typically grown in one of 
two ways: epitaxially on a sapphire substrate at high 
temperature, for a high-quality device on a limited scale, 
or, for large-area or flexible devices, on glass, plastic or 
metal, necessitating a sacrifice in quality. Therefore a 
fiexible, heat tolerant substrate is desired. In one de- 
sign, a GaN thin film is grown on graphene (functioning 
as the substrate and electrode) with the assistance of 
ZnO nanowalls to overcome the absence of GaN nucle- 



ation on pristine graphene ( [Chung et al] |2Q1Q[ ) . These 



graphene/GaN LEDs can then be transferred mechani- 
cally to other substrates. 

In the other devices, graphene is deposited onto pre- 
formed GaN films. The schematics and electrolumi- 
nescence intensity of two LEDs using In^jGai-a^N/GaN 
multi-quantum-wells (MQWs) as the emissive material 
are compared in Fig. [62^ and b. For the MLG anode 
device, the output power was shown to be comparable 
to an ITO-electrode device below 10 mA input current, 
becoming less efficient at higher input power. This is in 
part due to a larger forward voltage drop resulting from 
higher contact and sheet resistance of MLG compared to 
ITO. One way of correcting for high contact resistance 
is through doping. By intercalation of MLG semicon- 
ductor interfaces with bromine, it has been shown that 
tuning of the Schottky barrier height is possible over a 
wide range (jTongay et a/., |2011|. The Schottky bar- 



rier height is the rectifying barrier for electrical conduc- 
tion across a metal-semiconductor interface (here using 
graphene as the 'metal'). Bromine intercalation serves to 
increase interlayer separation between graphene sheets, 
as well as increasing free hole carrier density and decreas- 
ing the Fermi level due to electron transfer from carbon 
to bromine. 

Organic LEDs use polymers or small molecules for 
the emissive material. For efficient injection and charge 
conversion, the electrode work function must be well 
matched to the HOMO and LUMO of the emissive ma- 
terial. When coated with a thin and planarizing layer 
of PEDOT:PSS, ITO has an appropriate work function 
for hole injection as well as high conductivity and trans- 
parency; the latter is important so that the light gen- 
erated within the active layer can be emitted from the 
device. By blending an electroluminescent polymer with 
an electrolyte solution, (Matyba et a/., 2010) have con- 
structed an OLED device referred to as a light-emitting 
electrochemical cell (LEG). It uses entirely solution- 
processable carbon-based materials, eliminating the need 
for a metal electrode. The composition uses chemi- 
cally derived graphene for the transparent cathode, an 
active layer solution, and the conducting polymer PE- 
DOT:PSS as the transparent anode. The active layer so- 
lution consists of a blend of a poly (par a-phenylene viny- 



lene) copolymer called "Super Yellow" (SY) and an elec- 
trolyte in the form of the salt KGF3SO3 dissolved in 
poly (ethylene oxide) (PEO). The design of this device 
is shown in Fig. [62j 3. Because of the interfacial layer 
formed by the electrolyte, it mitigates the need for work 
function matching. In addition, the stability of the elec- 
trodes removes the possibility of ions from the electrode 
"leaking" into the active solution. 
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FIG. 62 Design and function of graphene LEDs and 

LECs. (a) Device schematics for inorganic LEDs incorporat- 
ing graphene as a growth substrate and cathode on the left, 
and as an anode on the right. With graphene forming the sur- 
face in contact with the substrate, it becomes easy t o move 
the entire device to a new substrate. Adapted from (Chung 



et al. 2010) and (Jo et al. 2010a). (b) Power- dependent elec- 
troluminescence of each device, showing strong blue emission 
spectra over a range of input currents. The left-hand device 
has been transferred to a plastic substrate, (c ) Schematic of 
an LEG. Adapted from |Matyba et al.\ |2010t . Blending of 
the emissive material with an electrolyte solution allows the 
rearrangement of charge at the electrode interface, facilitat- 
ing balanced electron-hole injection. As both electrodes are 
transparent, light is emitted through both surfaces. 



3. Photovoltaics 

Photovoltaic devices (PVs, also known as solar cells) 
are designed to harvest solar energy by producing cur- 
rent from materials that exhibit the photovoltaic effect. 
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Upon absorption of incoming light by the active material, 
electrodes collect the photogenerated charge carriers. Be- 
cause the electrodes physically surround the active mate- 
rial, at least one of them is transparent to allow light to 
reach the active layer. In some cases, both electrodes are 
transparent, though typically a metal such as aluminum 
or silver is used. The overall energy conversion efficiency 
T] of the PV is determined by the product of individual 
efficiencies: reflectance, thermodynamic, charge carrier 
separation and conductive. In practice, these may be dif- 
ficult to measure, so the overall efficiency is represented 
as: 



rj = 



p. 



Vqc X Isc X FF 
p. 



(33) 



where Isc is the maximum short circuit current, Vqc 
is the maximum open-circuit voltage and FF is the fill 
factor, defined as the ratio of actual maximum obtain- 
able power, given as Pmax = Vmax x Imax, to the theo- 
retical maximum, equal to Vqc x Isc- Pine is the total 
power of the incident light falling on the detector, and so 
77 represents the fraction of absorbed photons converted 
to current, and is also known as internal photocurrent 
efficiency. The best widely- available PVs are made from 
silicon, though comparable semiconductor PVs are in- 
creasingly available, such as thin-film cadmium telluride, 
copper indium selenide, and single-junction gallium ar- 
senide. Silicon solar cells have achieved an 77 value as 
high as 25%. The FF for commercial cells can reach 
70%, with 40% to 70% being typical. 

One instinctive direction is to integrate graphene with 
silicon wafers. By depositing a graphene sheet onto n- 
silicon, the authors of reference ( |Li et a/.||2010a| have cre- 
ated a Schottky junction, which provides faster switching 
speeds and better efficiency due to a lower forward volt- 
age drop. The graphene functions both as transparent 
upper electrode and antireflection coating, and through 
intrinsic electric field properties, aids in charge separa- 
tion and hole transport. Though the work is preliminary, 
opportunities for optimization are apparent, such as pas- 
sivation of the silicon to improve the interface - typically 
the deposition of an oxide layer forms a protective surface 
to prevent corrosion and reduce reactivity. 

One way of increasing the efficiency of devices contain- 
ing graphene films is to better match the work function 
of the graphene to the neighboring components, and sev- 
eral recent efforts have been made toward engineering 
the graphene work function (pc In designing a simple 
FLG/Si heteroj unction device, Ihm et al. revealed a di- 
rect relationship between the number of layers in the FLG 
and the resulting Vbc, with the best matching at 4 lay- 



ers, despite disordered stacking (Ihm et al. , 2010). As the 



Si substrate is unchanged, the Vqc is determined by the 
work function of the CVD-grown graphene layers, which 
was shown to better match the silicon as the number of 



layers increased, approaching the silicon work function 
(j)Si from below. 

Organic cells made from polymers can be manufac- 
tured at a lower expense, though they have not yet 
achieved comparable efficiencies to silicon. A popular 
choice of active layer for proof-of-concept devices is an 
organic polymer dispersion consisting of a narrow-band 
electron donor, poly(3-hexylthiophene) (P3HT), and a 
fullerene derivative, [6,6]-phenyl-C6i-butyric acid methyl 
ester (PCBM), which functions as an electron acceptor. 
Cells with this kind of donor /acceptor blend constitut- 
ing the active layer are referred to as bulk heteroj unction 
(BHJ) cells. As we will discuss, graphene has been inte- 
grated into inorganic and BHJ cells in various roles, as 
summarized in Table |V[ Most of the work on graphene 
PVs is at an early stage and in general they do not yet 
contend with the best performance obtained from simi- 
lar organic solar cells made with ITO electrodes. Most 
of the graphene PVs do not reach an r] of 2%, whereas 
those with ITO manage upwards of 6%. 



The recent work presented by (Wang et al., 2009a) 



illustrates some of the design considerations and opti- 
mization problems associated with a simple graphene- 
modified BHJ cell. CVD-grown FLG films function 
as the anode in their device, filling the role of ITO 
(or other TCE). The films varied in thickness from 
6-30 nm, with Rs from 1350 to 210 Vt/U and T 
from 91% to 72%. The complete device structure 
is: FLG graphene/PED0T:PSS/P3HT:PCBM/LiF/Al, 
and is shown in Fig. [63^ . The device is mounted on 
glass, which provides protection and structural integrity 
and permits light, which passes through the FLG layer 
to be absorbed by the P3HT:PCBM active layer, gener- 
ating the charge that will be collected by the FLG and 
LiF/Al electrodes. With a simple substitution of unmod- 
ified FLG film for ITO, the device efficiency was reduced 
by approximately 30- fold (from rj of 3.10% to 0.21%). 

This efficiency reduction was attributed to several fac- 
tors. First, the hydrophobicity of graphene precludes uni- 
form coating of the PEDOT:PSS layer, which is essential 
as a planarizing buffer layer for reducing surface rough- 
ness, as well the hole injection barrier between P3HT and 
graphene. To address this, the graphene was modified by 
UV/ozone treatment, which improved wettability by in- 
troduction of OH and C=0 groups to the surface. This 
resulted in improved Isc and Vqc values, with slightly 
reduced FF for an overall improvement to r] of 0.74%. 
The decrease in FF was attributed to increased series 
resistance of the cell, as disruption of the aromatic struc- 
ture by surface groups results in decreased conductivity. 

To circumvent this effect, the graphene was modi- 
fied instead by noncovalent functionalization with pyrene 
buanoic acid succidymidyl ester (PBASE). This resulted 
in the same improvement in Vqc and further improve- 
ment in both Isc and FF, resulting in overall 77 of 1.71%, 
reaching slightly more than half of the ITO version of 
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the device. The effect of PBASE modification was mul- 
tifaceted - in addition to avoiding 7r-conjugation disrup- 
tion and maintaining conductivity, the work function of 
graphene, ^c, was increased from 4.2 eV to 4.7 eV, in- 
creasing the open circuit potential as well as resulting in 
better matching between the Fermi level of FLG and the 
HOMO of PEDOT for more efficient hole collection. 
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FIG. 63 Schematics for some photovoltaic graphene 
devices, (a) Schematic and energy diagram of a bulk het- 
eroj unction cell that uses grap hene as a tra nsparent anode. 
Adapted from reference (Wang et al. 2009a|. (b) Construc- 



tion of an rGO-wrapped PDI hybrid wire, showing the in- 
teraction between graphene and PDI. Adapted from refer- 
ence (Wang et al, 2010). (c) Layered graphene/CdS, using 
graphene as the electron acceptor. By stacking several bilay- 
ers, the Isc and 77 of th e device are subs tantially improved. 
Adapted from reference (Guo et al. 2010). 



Many improvements to graphene BHJ cell designs have 
appeared recently, including using metals for doping or 
blocking. By immersing CVD-grown graphene films into 
AuCls, Shi et al. showed formation of Au particles on 
the graphene surface by spontaneous reduction of metal 



ions. The extent of the metal doping, controlled by im- 
mersion time, effectively tunes the surface potential, over 
a range of -0.5 eV (|Shi et aL\ |2Q10D. Ng et al. ex- 



plored the effect of introducing graphene to Ti02 nanon- 
structed films, using reduced graphene oxide (rGO, see 
Sec. IV. H) scaffolds to harvest photogenerated charge. 



and report up to a 90% enhancement in photocurrent 



over pristine Ti02 under UV illumination ( |Ng et al. 



2010). While pristine Ti02 films yield a maximum rj 



of 7.4%, the rG0-Ti02 nanocomposite films yield max- 
imum 77 of 13.9% for in situ rGO (photo-catalytic re- 
duction) and 11.4% for hydrazine rGO (chemical reduc- 
tion using hydrazine). In these devices, graphene serves 
as an efficient collector of generated electrons, and pre- 
vents recombination. It has also been shown that the 
introduction of a thin sol-gel processed titanium sub- 
oxide layer (TiO^c) for hole-blocking can improve de- 
vice efficiency 2-fold ( |Choe et al.\^oTo\) . In this CVD 
graphene/PEDOT:PSS/P3HT:PCBM/TiO^/Al device, 
the Isc and Vqc values closely match those for ITO, 
with only the lower FF preventing the overall efficiency 
from being comparable. 

An interfacial dipole layer (consisting of WPF-6-oxy-F, 
a synthesized polymer) in an inverted-structure organic 
solar cell can reduce the work function, increase built- 



in potential and improve charge extraction (Jo et al. 



2010b). CVD-grown graphene films were used as the 



device cathode for design flexibility and the ability to 
build stacked devices. The ionic or polar groups of the 
WPF-6-oxy-F form interface dipoles, better matching 
to the LUMO of PCBM (4.2 eV). Untreated MLG film 
has = 4.58 ± 0.08 eV, close to HOPG at 4.5 eV. It 
was reduced by 0.05±0.03 eV using poly(ethylene oxide) 
(PEO), 0.22±0.05 eV with CS2CO3, and 0.33±0.03 eV 
with WPF-6-oxy-F. 

Though graphene finds a natural function as a TCE, 
it also presents potential for additional utility by virtue 
of its chemical structure. In particular, its aromatic 
skeleton can function as a nucleating agent for the self- 
assembly of organic nanostructures. Using rGO flakes 
as a seeding agent for planar aromatic molecules (which 
tend to self-assemble), "hybrid" wires were formed 
through noncovalent hydrogen bonding and tt — tt stack- 
ing interactions ( Wang et al.\ 2010[ ). By choosing a 



photoluminescent molecule (in this case N,N'-dioctyl- 
3,4,9, 10-perylenedicarboximide, or PDI) which is effec- 
tively "wrapped" in graphene, a photovoltaic active 
material can be produced, as shown schematically in 
Fig. [63) 3. Hydrothermal reduction forms PDI-graphene 
hybrid wires (PDI-G) which can be integrated into a PV 
device, with properties as listed in Table [V| 



I 

4. Graphene quantum dots (GQDs) and graphene-QD 2D SLG. Quantum dots are described as semiconduc- 

nanocomposites 

Another interesting development is the formation of 
OD graphene nanoparticles or nanowells from isolated 
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Device construction 




Voc 


FF 




Role of 


Reference 




(mA/cm ) 


(V) 


(%) 


(%) 


graphene 




IT0/PED0T:PSS/P3HT:PCBM/LiF/Al 


9.03 


0.56 


61.1 


3.10 


Transparent anode ( 


Wang et a/.||2009a| 


Graphene/n-Si/Ti/Pd/Ag 


6.5 


0.48 


55 


1.65 Transparent cathode 


jLi et a/.l|2010al 


CVD graphene/PEDOT/CuPc/Ceo/BCP/Al 


4.73 


0.48 


52 


1.18 


Transparent anode iDe Arco et ah, 201C 




IT0/PED0T:PSS/P3HT:graphene/LiF/Al 


4.0 


0.72 


38 


1.1 


e acceptor 


pu et a/., 2009) 


IT0/PED0T:PSS/P3HT:PCBM-graphene/Al 


5.3 


0.64 


41 


1.4 


Dopant 


\u\i et a/., 2010) 


Graphene/PED0T:PSS/P3HT:PCBM/LiF/Al 


6.05 


0.55 


51.3 


1.71 


Transparent anode (Wang et a/., 2009a 


L 


MLG/PEDOT:PSS/P3HT:PCBM/TiOa./Al 


9.03 


0.60 


48 


2.60 


Transparent anode 


jChoe et al. 2010^ 


Al/PEDOT:PSS/P3HT:PCBM/WPF-6-oxy-F/MLG 


6.61 


0.57 


33 


1.23 Transparent cathode 


fJo et a/., 2010^ 


IT0/PED0T:PSS/P3HT:GQDs/Al 


6.33 


0.67 


30 


1.28 


e~ acceptor 


(Li et a/., 2011) 


ITO /PEDOT /P3HT:PDI-G /LiF / Al 


3.85 


0.78 


35 


1.04 


e~ acceptor 


(Wang et a/., 2010) 




Layered graphene/CdS QDs on ITO glass* 


1.08 


0.68 




16 


e~ acceptor 


(|Guo et a/.||2010[ 



TABLE V Function of graphene photovoltaic devices. A summary of some recently developed photovoltaic devices that 
have incorporated graphene in various capacities, compared to an ITO device (shown at the top). 
*For 8 bilayers. 



tors (typically nanoparticles) with excitons confined in 
3 dimensions, which results in interesting properties that 
fall somewhere between bulk semiconductors and discrete 
molecules. Most inorganic QDs (CdSe, CdTe, InP) are 
fluorescent, with a large Stokes shift and narrow, size- 
dependent emission peaks, making them useful in fluo- 
rescence applications. 

Recently, a straightforward electrochemical method 
of producing GQDs from freestanding graphene film 
was reported (Li et al. 2011). By treating with O2 



and using the graphene as the working electrode in 
a cyclic voltammetry cell, they were able to produce 
hydrophilic, monodisperse GQDs 3-5 nm in diameter. 
Interestingly, they show an excitation-dependent emis- 
sion wavelength and intensity, distinct from other QD 
compositions. With high-resolution XPS, it was deter- 
mined that the as-prepared particles' surface contained 
hydroxyl, carbonyl and carboxylic acid groups, both 
making them water-soluble, as well as paving the way 
for potential conjugation to other molecules. They are 
also soluble in some organic solvents for use in the non- 
aqueous phase of devices. With this in mind, they were 
integrated into polymer PV cells, with a composition 
of IT0/PED0T:PSS/P3HT:GQDs/Al to form a metal- 
insulator- metal (MIM) device, giving r] = 1.28% after 
some optimization, outperforming the GQD-free assem- 
bly. 

Graphene sheets have also been integrated into QD 
nanocomposites. Such a device been developed with CdS 
QDs and graphene, harnessing the electrical conductiv- 



ity of graphene to form an electron-transport matrix ( Cao 
\et a/.[|2QlQ| ). The fabrication process uses graphene oxide 
in dimethyl sulfoxide (DMSO), allowing GO reduction 
and CdS deposition simultaneously. They report picosec- 
ond electron transfer from the excited QDs to graphene 
using time-resolved fluorescence microscopy. Another 




(b) 




300 400 500 

Wavelength (nm) 



350 400 450 500 550 6OO 650 700 750 

Wavelength (nm) 



FIG. 64 Absorption and photoluminescence of 
graphene quantum dots in water, (a) Absorbance spec- 
trum of GQDs, showing a broad UV absorption. Inset shows 
a vial of GQDs in water, (b) Excitation-dependent emission 
of GQDs. As the excitation wavelength is red-shifted, the 
emission is also red-shifted to a lesser extent, along with di- 
minished intensity. 



promising PV device consists of layers of graphene and 



CdS QDs built upon an ITO glass support (Guo et al. 



2010), as shown in Fig. 63 3. The effect of the number of 



graphene/QD bilayers on the photocurrent and voltage 
was investigated, and a substantial increase in the Isc 
up to 8 layers was shown, while the Voc held around 0.68 
V for 2 or more layers. An rj value of 16% is reported, 
considerably surpassing previous versions of carbon/QD 
solar cells {rj < 5%). This is attributed to the thin- 
layered structure being effective at charge collection and 
transport. 



XI. GRAPHENE SENSORS 

In addition to the optoelectronic applications discussed 
above, important research has also been conducted on 
graphene sensors. After the isolation of graphene in 2004 



(Novoselov et al. 2004), graphene-based nanomaterials 
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have for the most part replaced CNTs in sensor apphca- 



tions (Geim and Novoselov 2007 Li and Kaner 2008), 



due to higher sensitivity to adsorbed materials (Liang 



and Zhi 



2009 



Pumera 



2009) 



et all '2009; "Fowler et al 


|2009| ILu et aL\ |2009t Qazi 


et al] 2007jjRobinson et al 


2008a jSchedin et al , 2007b). 


The authors of ( 


Robinson et a/., 


2008a 


) demonstrated a 



A. Electrochemical sensors 

Graphene exhibits improved electrochemical response 
when compared with other electrodes such as glassy car- 



molecular sensor based on rGO capable of detecting toxic 
gases with ppb sensitivity. In comparison to CNT based 
sensors, RGO based sensors show similar performance 
with greatly reduced noise. The detection of NO2 as 
low as 60 ppb as been demonstrated by measuring the 



changes in conductivity of flakes in thin graphite (Qazi 



2009), and CNTs (Alwarappan et al 



bon electrodes (Shang et al. , 2008), graphite (Shan et al. 



2009b). The authors of (Zhou et al. 



2009 Wang et al. 



et al. 2007). Detection of organic vapors like nonanol, 
octanic acid, and trimethyl amine was reported in ref- 
erence (Dan et a/., 2009). Furthermore, the authors of 



that graphene exhibits a wide electrochemical potential 
window of ca. 2.5 V in 0.1 M phosphate-buffered saline 



2009) have shown reference (ISchedin et al. 



2007b ) achieved single molecule 



solution (PBS) at pH 7.0. In Fig. 65, graphene platelets 
(GNPs) were shown to simultaneously detect the small 
organic molecules dopamine, ascorbic acid, and uric acid. 
Other groups have explored the electrochemical proper- 
ties of reduced graphene oxide (rGO) and graphene-based 



nanomaterials (Wang et al.\ 2009bl IZhou et a/., 2009). 



These studies used different probes such as nucleic acids, 
potassium ferricyanide, dopamine, and acetaminophen. 




FIG. 65 Graphene platelets (GNPs) for electrochemi- 
cal sensing (A) (a) and (b) Transmission electron microscopy 
(TEM) images of GNPs at different magnifications; (c) high- 
resolution TEM image of GNPs showing a nanoflake with a 
knife-edge or conical structure with open graphitic planes; (d) 
energy-dispersive X-ray spectrum, showing the chemical com- 
position of GNP films. (B) (a) and (b) Cyclic voltammetry 
profiles of GNPs and bare GCEs, respectively, in a solution of 
50 mM PBS buffer at pH 7, containing 1 mM ascorbic acid , 
0.1 mM dopamine, and 0.1 mM uric acid ( |Shang et al. 2008). 



CNTs have shown the ability to detect gases when 



decorated with nanoparticles (Kong et al. 2000 Sippel- 



Oakley et al. , 2005 


Star et al. 




2004). Graphene based 


sensors also exhibit an electrica 


. response to gaseous CO, 



H2O, NO2, and NH3. Sensors based on graphene oper- 
ate by measuring a change in resistivity resulting from 



detection with an electrical sensor made of few-layered 
graphene. 





FIG. 66 Electrical response of graphene to various 
gases (a) Concentration An of chemically induced charge car- 
riers in single-layer graphene exposed to different concentra- 
tions C of NO2. Upper inset: SEM micrograph of this device; 
the width of the Hall bar is 1 /xm. Lower inset: the changes 
in resistivity p and Hall resistivity pxy of the device with gate 
voltage Vg. The ambipolar field effect is clearly illustrated, 
(b) Changes in resistivity at zero magnetic field as caused by 
exposure of graphene to different gases diluted in a carrier 
gas of either He or N2 (final concentration 1 ppm). Region I: 
device is in vacuum; II: analyte exp osure; HI: analyte evac u- 
ation; and IV: annealing at 1500°C ([Schedin et"arj |2007b ). 



Fig. 66 shows the electrical response of graphene to 
various gases. Here the donation of electrons to graphene 
increases the resistance of the device, whereas scaveng- 
ing of electrons decreases the overall resistance. In this 
case, Schedin and co-workers managed single molecule 
detection with pure graphene due to an unintentional 
functionalization of their device with a residual polymer 
layer from lithographic resist. The response of the sensor 
was measured before and after the removal of polymer 
resist flSchedin et al.\ |2007b|). Without the resist, the 



sensitivity of the device dropped by one to two orders 
of magnitude (Dan et al. 2009) . Impurities and defects 



were also shown to affect the response to gas adsorption 



the adsorption of gas molecules (Arsat et al. , 2009 Dan onto graphene (Ao et al. 2008 Zhang et al. 2009b). 
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B. Graphene as a biosensor 

The same mechanism ahows graphene and its deriva- 
tives to be used in biosensors. The electrochemical re- 
sponse of hydrogen peroxide (H2O2) on an rGO elec- 
trode has been studied which resulted in improved per- 
formance compared to glassy carbon, graphite and CNTs 
(Zhou et al, 2009). A similar effect was observed in 



the electrochemical behavior of the small molecule reduc- 
ing agent nicotinamide adenine dinucleotide (NADH) on 
graphene- modified electrodes (Tang et al. |2009p . In both 
the cases, the superior electrolytic activity of graphene 
is attributed to the high density of edge-plane-like de- 
fective sites on graphene. The enhanced oxidation of 
NADH on graphene-modified electrodes is strongly con- 
firmed when compared with bare edge plane pyrolytic 
graphite electrode (EPPGE). High density of the edge- 
plane-like defective sites contributes to enhanced oxi- 



dation/reduction of biomolecules (Banks et al. 2004). 
A multilayer graphene nanofiake film electrode demon- 
strated highly resolved simultaneous detection of uric 
acid, ascorbic acid and dopamine with the detection 
limit as low as 0.17 /iM (Zhou et al. 2009). It was re- 



ported that graphene exhibited better sensing capabil- 
ity of dopamine than CNTs in resolving ascorbic acid, 
dopamine and serotonin (Alwarappan et al. 2009|). In 
a similar study, graphene exhibited high sensitivity to 
dopamine with a linear range of 5-200 /iM and a bet- 
ter performance compared to multi-walled CNTs ( |Wang 
et al., 2009b). This performance is attributed to high 
conductivity, large surface area and tt — tt bond interac- 
tion between dopamine and graphene. 

An electrochemical sensor based on functionalized 
graphene as been designed for the sensitive detection of 
acetominophen (Kang et al., 2010). The electrochemical 



behavior of acetominophen on graphene modified GCEs 
was studied with the help of cyclic voltammetry and 
square wave voltammetry. Chen et al. designed a glucose 
biosensor using graphite platelets based on a composite 
material including GNPs, Nafion binder, and glucose ox- 
idase (Fu et al. 2009j). GO has been successfully shown 



to support electrical writing of the redox centers of vari- 
ous metalloproteins, without altering their structural in- 



tegrity or biological activity (Fu et al. 2009). 

Biosensors produced from GO and from graphene 
amine made by treating GO with nitrogenous plasma or 



ethylenediamine were developed in reference (Mohanty 
and Berry, 20081). Once treated, graphene layers were 
electrostatically adsorbed on to a silica substrate. Elec- 
trical measurements showed that rGO sheets acted as 
p-type semiconductors with high resistance and with 
extremely low carrier mobilities. Fig. [67| shows nega- 
tively charged bacteria adsorbed on positively charged 
graphene amine. This adsorption of single bacteria is ac- 
companied by a 42% increase in conductivity. In another 
experiment, single- stranded DNA (ssDNA) was chemi- 



cally grafted onto a graphene surface. Hybridization of 
the tethered strands by fiuorescently-tagged complemen- 
tary strands was observed. The authors further observed 
that the tethering of ssDNA on graphene was preferred 
on thicker sheets and on wrinkled surfaces of rGO rather 
than on fiat surfaces. Initial tethering of ssDNA more 
than doubled the composite's conductivity, implying an 
increase in hole density at the graphene layer on the or- 
der of 5.61 X 10^^ cm~^. Monolayers of ssDNA molecules 
were adsorbed on both sides of graphene sheets by tt — tt 



stacking (Patil et al. 2009). The experimental procedure 



used to produce ssDNA coated graphene sheets involved 
vigorous sonication of a suspension of ssDNA and hy- 
drophilized graphene. 




G.3iDfJ A ^Chdiybridlzcd) 
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FIG. 67 Negatively charged bacteria adsorbed on pos- 
itively charged graphene amine (a) Electrical response of 
graphene amine (GA) device upon attachment of a single bac- 
terium on the surface of GA (inset 1). LIVE/DEAD confocal 
microscopy test on the bacteria deposited on GA confirmed 
that most of the bacteria were alive after the electrostatic 
deposition (inset 3). a) Alive and d) Dead, (b) DNA transis- 
tor: ssDNA tethering on GO increases the conductivity of the 
device. Successive hybridization and dehybridization of DNA 
on the G-DNA device results in completely reversible increase 
and restoration of conductivity. Inset shows a G-DNA(ds) 
sheet with w rinkles and folds clearly visible (Mohanty and 
Berry||2008|. 



Graphene nanomaterials have also been incorporated 
into functionalized biosy stems integrated with DNA, 
peptides, proteins, aptamers, and cells. Ohno et al. re- 
ported the use of a graphene-electrolyte field effect tran- 
sistor to detect dissolved bovine serum albumin (BSA), 
which is an important protein for biochemical applica- 



tions ( |Ohno et al. , 2009). The sensor was made of me- 
chanically exfoliated single layer graphene supported on 
a silicon dioxide wafer. An electrolyte solution and ref- 
erence electrode completed the circuit. The authors ob- 
served the electrical response of their device with ad- 
sorption of BSA from standard solutions. The detection 
limit of BSA was as low as 0.3 nM. The same research 
group fabricated an apt amer- modified graphene-FET im- 
munosensor. Functionalization of the G-FET with ap- 



tamers was confirmed by atomic force microscope ( Ohno 



et 



2010|. 
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XII. CONCLUSION 

Graphene's journey is quite remarkable considering the 
tremendous scientific and technological impact this ma- 
terial has had on the scientific community. From the 
first study of graphene's band structure by Wallace in 
1947 at McGill University, in the footsteps of whom the 
authors of this review have followed, to the pioneering 
experimental studies that were awarded the Nobel Prize 



in Physics in 2010. The former ( |Wallace[ |1947| ) led to 
the theoretical foundation of graphene, the unique and 
simple dispersion described by the Dirac Hamiltonian, 
while the latter work culminated with the more recent 
isolation of single graphene crystals in 2004. This last 
discovery ignited an explosion of theoretical and experi- 
mental work on graphene and inspired the writing of this 
review of graphene's remarkable experimental properties 
with a practical emphasis. The decision was made to fo- 
cus on the experimental properties, since we believe that 
they make for an excellent introduction on many different 
experimental techniques which are also relevant to other 
fields in condensed matter physics, materials engineering, 
chemistry and biophysics. 

Writing a review is necessarily selective and the goal 
was not to be fully comprehensive, but to include the 
topics the authors believe to be the most interesting for 
the widest possible range of new researchers in this field. 
This naturally includes particular highlights such as the 
ambipolar effect for transistors, the anomalous quantum 
Hall effect, which brought this field to room tempera- 
ture, the remarkable optical properties, while remaining 
a good conductor, the fascinating mechanical properties 
and selective sensitivity to the environment, with great 
outlook as sensor material. All of this would be use- 
less without the possibility to synthesize and fabricate 
graphene, which is extensively reviewed here. 

History has shown that predicting the future of the 
field and all its possible applications is a difficult and 
risky endeavor. As the experience of the authors is very 
modest in that respect, we shall only say that a lot re- 
mains to be done and that graphene will be an important 
topic for many years to come. 
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